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ABSTRACT  OF  THE  THESIS 


The  Effect  of  Methanol  on  the 
Mobilization  of  Tetrachloroethylene 
in  a  Water-Saturated  Soil  Column 
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The  contamination  of  surface  waters  and  groundwater  by  chlorinated  solvents 
continues  to  be  a  major  environmental  issue.  Today,  the  challenge  to  environmental 
engineers  is  not  necessarily  which  remediation  technique  to  use,  but  how  to  enhance  the 
remediation  technique  in  order  to  accelerate  the  clean-up  process. 

The  "pump-and-treat"  method  for  remediating  organic  contaminants  is 
experimenting  with  the  use  of  chemical  additives  to  increase  the  efficiency  rate  of 
subsurface  contaminant  removal.  Because  of  their  favorable  chemical  properties,  alcohols 
such  as  methanol  are  often  used  in  conjunction  with  water  to  flush  a  contaminated  area, 
and  subsequently  pump  it  to  the  surface  for  treatment.  For  this  study. 
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the  objectives  were:  (1)  to  reconfinn  the  impact  of  a  cosolvent,  such  as  methanol,  on  the 
observed  retardation  of  nonpolar  organic  solute  (tetrachloroethylene,  PCE)  transport  in  a 
water-saturated  porous  medium;  and,  (2)  to  provide  evidence  that  ineffective  mixing  of 
two  inhomogeneous  fluids,  such  as  methanol  and  water,  can  negatively  impact  the  desired 
cosolvent  effect. 

For  the  one-dimensional  laboratory  column  at  the  flowrate  tested  (0.745  mL/min), 
methanol  failed  to  produce  the  hypothesized  cosolvent  effect.  The  degree  of  retardation 
of  PCE  did  not  decrease  log-linearly,  but  instead  showed  no  significant  difference  with 
increasing  methanol  fraction.  Therefore,  the  column  system,  under  the  conditions  used  in 
this  study,  may  have  been  subject  to  the  ineffective  fluid  mixing  potentially  associated 
with  cosolvent-water  systems.  However,  breakthrough  behavior  for  methanol  was 
generally  symmetric  and  failed  to  exhibit  any  signs  of  ineffective  mixing.  The  nature  of 
these  results  vmderscores  the  complexity  of  the  problem  of  mixing  of  inhomogeneous 
fluids. 
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Chapter  L  Introduction 

Chlorinated  solvents,  used  as  cleaning  agents  in  a  vast  array  of  industries,  from 
dry  cleaning  to  aerospace,  are  known  or  suspected  as  either  carcinogens  or  mutagens. 

The  frequent  industrial  use  of  these  chemicals,  such  as  tetrachloroethylene  (PCE), 
trichloroethane  (TCA),  and  trichloroethylene  (TCE),  has  led  to  widespread  contamination 
of  soil  and  groundwater.  Because  of  the  high  volatility  of  these  organic  compounds,  their 
concentrations  in  surface  waters,  which  are  only  a  few  micrograms  per  liter,  is  not  a  water 
quality  issue.  However,  it  is  in  the  groundwater  where  the  concentrations  of  chlorinated 
solvents  will  be  thousands  of  times  higher  since  the  contaminants  cannot  volatilize,  and 
hence  will  tend  to  remain  in  water  or  sorb  to  soil.  Most  frequently,  the  corrective  action 
for  subsurface  contamination  is  to  flush  the  contaminated  area  with  water  and 
subsequently  pump  it  to  the  surface  for  treatment.  This  conventional  "pump-and-treat" 
(Mackay  and  Cherry,  1989)  strategy  is  effective  for  removal  of  a  majority  of  the 
contaminant.  Unfortunately,  the  remaining  contaminant  is  still,  most  often,  well  above 
prescribed  limits. 

The  remaining  contamination  is  often  due  to  nonaqueous  phase  liquids  (NAPLs) 
in  the  subsurface.  NAPLs  are  immiscible  in  water  and  can  have  densities  that  are  either 
greater  than  water  or  less  than  water.  When  a  NAPL  has  a  density  greater  than  water,  it  is 
commonly  referred  to  as  a  dense  nonaqueous  phase  liquid  (DNAPL).  Conversely,  if  the 
NAPL's  density  is  less  than  water,  it  is  commonly  referred  to  as  a  light  nonaqueous  phase 
liquid  (LNAPL).  Examples  of  DNAPLs  include  chlorinated  solvents  such  as 
trichloroethylene.  Some  examples  of  LNAPLs  include  gasoline,  petroleum  oil,  and  diesel 
fuel.  The  continued  presence  of  NAPLs  in  the  subsurface  poses  a  major  clean-up 
problem.  It  is  critical  to  locate  and  remove  any  NAPLs  from  a  hazardous  waste  site 
because  NAPLs  provide  a  long-term  source  of  pollution  due  to  their  low  solubilities. 
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Because  of  the  limitations  of  the  "pump-and-treat"  method  for  remediating  NAPL  organic 
contaminants,  alternative  remediation  techniques  such  as  bioremediation  and  chemical 
additives  are  being  considered. 

1.1  Bioremediation 

Bioremediation  is  the  use  of  microorganisms  to  degrade  subsurface  contaminants. 
The  microorganisms,  generally  bacteria,  convert  harmful  chemical  compounds  to  less 
harmful  chemical  compounds  in  order  to  effect  remediation  of  a  contaminated  site.  To 
survive,  the  microbes  require:  a  carbon  supply,  an  energy  supply,  and  nutrients.  The 
types  and  amounts  of  carbon,  energy,  and  nutrients,  are  part  of  the  overall  design  of  the 
remediation  system.  Therefore,  a  detailed  understanding  and  control  of  the  site 
hydrogeology  is  required  to  implement  an  efficient  remediation  effort.  Contaminant 
factors,  such  as:  solubility,  volatility,  viscosity,  and  toxicity;  and  soil  conditions,  such  as: 
permeability,  soil  type,  depth  to  grovmdwater,  mineral  content,  oxidation/reduction 
potential,  and  pH;  may  have  a  profound  affect  on  the  general  remediation  design. 

The  key  to  successful  bioremediation  is  engineering  the  system  such  that  all  the 
microbe  requirements  can  combine  in  the  subsurface.  However,  engineering  the  system 
is  not  as  simple  as  one  might  expect.  The  geologic  medium  for  remediation  must  be 
permeable  enough  to  allow  the  introduction  of  energy  and  nutrients.  Many  intraparticle 
pores  are  too  narrow  to  allow  microbe  access.  If  the  microbes  cannot  reach  the  zone  of 
contamination,  the  remediation  effort  will  fail.  Another  limitation  to  bioremediation  is 
that  it  is  difficult  to  biostimulate  immobile  zones.  Most  often,  the  contaminant  must 
diffuse  out  to  the  biostimulated  zone  before  bioremediation  can  occur.  Finally, 
bioremediation  is  still  in  the  developmental  stage  and  there  are  few  documented  cases  of 
using  bioremediation  in  the  saturated  zone. 
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1.2  Chemical  Additives 

Enhanced  remediation  techniques  using  chemical  additives  are  being  introduced 
to  increase  the  efficiency  rate  of  contaminant  removal  from  the  subsurface  (Palmer  et  al, 
1992;  Augustijn  et  al,  1994).  In-situ  solvent  flushing  involves  injecting  a  mixed  solvent 
into  the  subsmface  of  a  site  contaminated  with  organic  chemicals.  Surfactants,  alcohols, 
or  mixtures  of  these  chemicals  are  often  used  as  the  solvent  (Imhoff,  1995).  When  the 
solvent  is  in  sufficient  quantities  relative  to  water,  it  is  often  referred  to  as  a  water- 
miscible  organic  cosolvent.  Miscibility  describes  a  compound's  ability  to  mix  in  any  ratio 
with  water  without  separation  into  two  phases,  while  immiscibility  means  a  compound  is 
unable  to  mix  with  water.  Because  of  their  favorable  chemical  properties,  alcohols,  such 
as  methanol,  have  been  investigated  as  potential  cosolvents. 

1.3  Thesis  Overview 

The  objective  of  this  thesis  study  is  to  test  the  ability  of  methanol  as  an  effective 
cosolvent.  Through  a  series  of  experiments,  this  research  will  investigate:  the  transport 
and  dissolution  of  methanol  in  a  water-saturated  porous  medium  and,  the  effects  of 
methanol  on  the  mobilization  of  PCE.  The  hypothesis  of  this  study  is  that  methanol  will 
decrease  the  retardation  factor  of  PCE  in  a  porous  medium,  thereby  confirming  the 
impact  of  cosolvents  on  remediating  hazardous  organic  compounds. 

In  Chapter  Two,  background  material  relating  to  organic  contaminants  will  be 
presented.  Chapter  Three  will  delve  into  the  theory  of  cosolvency  and  subsurface 
remediation.  In  Chapter  Four,  the  research  objectives  and  approach  for  this  thesis  will  be 
outlined.  Chapter  Five  will  detail  the  experimental  methods  used,  and  Chapter  Six  will 
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highlight  the  results.  Finally,  in  Chapter  Seven,  a  summary  of  the  findings,  and  the 
conclusions  reached  in  this  study,  will  be  presented. 
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Chapter  2.  Background 

Before  reviewing  the  underlying  theory  of  cosolvency,  it  is  instructive  to  examine 
background  material  relating  to  organic  contaminants  and  their  properties,  properties 
affecting  the  transport  of  a  NAPL-based  solute,  and  physical  properties  affecting  the  fate 
and  transport  of  organic  chemicals  in  groundwater. 

2.1  Organic  Contaminants  and  their  Properties 

Organic  contaminants,  such  as  TCE  and  PCE,  are  just  two  of  the  myriad  of 
chlorinated  solvents  found  in  the  subsurface.  The  extent  to  which  these  solvents 
contaminate  the  subsurface  depends,  in  part,  on  the  properties  of  the  organic 
contaminants.  Polarity  refers  to  the  extent  of  the  dipole  moment  in  a  molecule. 
Hydrophobicity  is  the  tendency  of  organic  compounds  dissolved  in  groundwater  to 
adsorb  onto  solid  surfaces.  Less  polar  molecules  are  more  hydrophobic.  Solubility  refers 
to  the  concentration  of  a  species  in  a  saturated  aqueous  solution.  TCE  and  PCE  have  low 
solubilities  due  to  their  lack  of  polarity;  hence,  they  are  commonly  referred  to  as 
hydrophobic  contaminants. 

2.2  Properties  Affecting  Transport  of  a  NAPL-based  Solute 

Several  key  properties  affect  the  transport  of  a  NAPL-based  solute.  Solubility 
plays  a  major  role.  Another  important  property  is  sorption.  Sorption  is  the  term  used  to 
describe  adsorption  and/or  absorption  when  independent  determination  of  the  processes  is 
impossible.  Adsorption  is  the  net  accumulation  of  a  solute  at  the  solid-liquid  interface; 
while  absorption  occurs  when  a  solute  is  transferred  from  one  phase  (i.e.,  water)  to 
another,  such  as  soil  or  organic  matter. 
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2.3  Physical  Processes  Affecting  Fate  and  Transport  of  Organic 
Chemicals  in  Groundwater 

Organic  contaminants  may  enter  the  subsurface  by  spills,  leaks,  or  intentional 
disposal.  In  the  saturated  zone,  these  organic  chemicals  will  either  dissolve  in  the 
groundwater  or  will  remain  immiscible  with  groundwater  in  the  form  of  NAPLs.  A 
review  of  the  physical  processes  that  most  strongly  influence  this  phenomena  will  be 
presented  to  assist  in  the  understanding  of  organic  contaminant  transport  and  distribution 
in  the  saturated  zone. 

2.3.1  The  Advection-Dispersion  Equation 

Transport  phenomena  in  the  saturated  zone  is  governed  by  the  Advection- 
Dispersion  equation  (e.g..  Freeze  and  Cherry,  1979): 


dc  =  a  (Dhd  5C)  -  a  (vx  ac)  (2-i) 

dt  dx  dx  dx  dx 

(1)  (2)  (3) 

where  Dhq  =  hydrodynamic  dispersion  coefficient  [L2/T] 

Vx  =  the  average  pore  water  velocity  [L/T] 
t  =  time  [T] 

X  =  distance  [L] 

C  =  solution  concentration  [M/L^] 

The  first  term  (1)  is  the  net  accumulation.  Term  (2)  represents  the  transport  due  to 
dispersion  and  term  (3)  represents  the  transport  due  to  advection. 

A  more  complete  modeling  of  the  transport  of  the  solute  in  the  aquifer  would 
consider  the  immobile  zones  in  each  molecule  of  porous  media.  To  accoimt  for  this 
property  of  the  media  the  concept  of  retardation  will  be  introduced  in  Section  2.3.4. 
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2.3.2  Advection 

Saturated-zone  transport  of  dissolved  contaminants  is  due  primarily  to  advection. 
Advection  is  the  process  by  which  a  solute  is  transported  by  the  movement  of  the  fluid 
itself  (i.e.,  groundwater  flow).  The  flow  of  groundwater  is  driven  according  to  Darcy's 
Law.  It  is  important  to  be  able  to  estimate  the  rate  at  which  groundwater  is  moving 
through  an  aquifer.  This  rate,  the  Darcy  velocity,  can  be  found  by  solving  for  the  velocity 
using  Darcy's  Law  (2-2). 

Q=-KA2h  (2-2) 

dL 

where  Q  =  flow  rate  [L^/T] 

K  =  hydraulic  conductivity  or  coefficient  of  permeability  [L/T] 

A  =  cross-sectional  area  [L^] 

0h  =  hydraulic  gradient 
dL 

The  Darcy  velocity  is  not  the  actual  groundwater  velocity  because  the  cross-sectional 
area,  given  by  A  in  (2-2),  of  an  aquifer  is  made  up  of  solids.  The  actual  cross-sectional 
area  is  much  smaller,  and  hence,  the  actual  groundwater  velocity  is  considerably  faster 
than  the  Darcy  velocity. 

Groundwater  will  flow  from  points  of  high  piezometric  head  to  points  of  low 
piezometric  head.  The  difference  in  water  elevation  and  pressure  head  between  points 
determines  the  hydraulic  gradient  which  drives  the  flow  of  groundwater.  In  gentle 
topography  and  under  unconfined  conditions,  the  gradient  normally  follows  the 
topography,  and  the  groundwater  flow  rate  can  be  estimated  vdthin  a  factor  of  10 
(Roberts  et  al,  1986).  However,  several  factors,  such  as  a  porous  media's  hydraulic 
conductivity,  can  make  the  determination  of  the  groundwater  flow  rate  quite  difficult  to 
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estimate.  Hydraulic  conductivity  is  a  measure  of  the  ability  of  a  porous  media  to  transmit 
a  fluid  (i.e.,  water)  and  is  dependent  on  the  properties  of  the  porous  media  and  the  fluid. 

A  subsurface  porous  media  found  in  nature  is  generally  heterogeneous-,  that  is,  the 
hydraulic  conductivity  varies  from  place  to  place.  The  porous  media  is  considered  to  be 
homogeneous  when  the  hydraulic  conductivity  is  the  same  throughout.  Sometimes 
hydraulic  conductivities  are  a  function  of  the  flow  direction.  Generally,  hydraulic 
conductivities  are  higher  in  the  horizontal  direction  than  in  the  vertical  direction  (Freeze 
and  Cherry,  1979).  Aquifers  that  have  the  same  hydraulic  conductivity  in  any  flow 
direction  are  called  isotropic,  while  those  aquifers  in  which  the  conductivities  varies 
according  to  direction  are  called  anisotropic.  Because  homogeneous  and  isotropic 
aquifers  are  rarely  found  in  nature,  one  must  collect  empirical  measurements  to  obtain 
reasonably  accurate  estimates  of  the  groundwater  velocity  in  a  flow  field. 

For  one-dimensional  soil  column  experiments,  such  as  the  one  used  for  this 
study,  a  homogeneous  and  isotropic  aquifer  can  be  simulated.  Because  a  pore  water 
velocity  is  typically  known  for  this  type  of  system,  the  flow  equation  (2-2)  is  not 
employed.  A  one-dimensional,  homogeneous  and  isotropic  porous  medium  will  limit 
variations  from  occurring.  However,  the  limitation  of  this  lab  system  is  readily  apparent 
in  that  it  is  not  representative  of  an  actual  environmental  system. 

2.3.3  Dispersion 

Dispersion  is  the  tendency  of  a  solute  to  spread  during  transport.  The  two  main 
causes  of  dispersion  are  mechanical  mixing  and  molecular  diffusion  (Freeze  and  Cherry, 
1979).  Mechanical  mixing  is  mixing  due  to  velocity  gradients.  The  factors  which  affect 
the  degree  of  mechanical  mixing  include;  different  sized  pore  channels  within  the  porous 
media,  variations  in  pore  geometry  (tortuosity),  and  fluctuations  in  the  local  flow 
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velocities  relative  to  the  mean  flow  direction.  Molecular  diffusion  is  the  net  flux  of  a 
solute  from  a  zone  of  high  concentration  to  a  zone  of  low  concentration.  Dispersion 
results  in  the  dilution  of  contaminant  pulses  and  the  attenuation  of  concentration  peaks. 
Therefore,  the  maximum  concentrations  diminish  with  increasing  distance  from  the 
source  (Mackay  et  al,  1985). 

Dispersion  can  occur  longitudinally,  in  the  principal  direction  of  flow;  or 
transversely,  perpendicular  to  the  principal  direction  of  flow.  The  hydrodynamic 
dispersion  coefficient  is  a  tensor  (Bear,  1972): 


D  =  aLDl+aTD2  +  D*  (2-3) 


where  D  =  hydrodynamic  dispersion  coefficient  [L^/T] 

D*  =  coefficient  of  molecular  diffusion,  which  is  generally  ignored  [L^/T] 
ttL  =  longitudinal  dispersivity  of  the  porous  medium  [L] 
aj  =  transverse  dispersivity  of  the  porous  medium  [L] 

For  the  one-dimensional  case. 


D*  =  V  and  D^  =  0  where  v  =  Darcy  velocity  [L/T] 


Therefore, 


D  =  ttL  (v) 


(2-4) 


Despite  the  practical  importance  of  the  dispersion  process,  there  is  currently  no 
method  to  confidently  predict  the  magnitude  of  dispersion  for  a  previously  unstudied 
field  situation.  For  simple  hydrogeological  systems,  the  spreading  is  believed  to  be 
proportional  to  the  flow  rate.  For  more  complex  systems,  the  dispersivity  appears  to 
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depend  on  the  structure  of  the  geologic  medium  such  that  it  varies  with  the  distance 
traversed  (Mackay  et  al,  1985). 

2.3.4  Retardation 

Retardation  is  the  apparent  slowing  of  contaminant  transport  relative  to 
groundwater  flow  due  to  sorption.  The  higher  the  fraction  of  the  contaminant  sorbed,  the 
more  retarded  is  its  transport  (Mackay  et  al,  1985). 

Sorption  of  halogenated  chemicals  to  soils  is  a  nonspecific  combination  of  two 
factors.  The  first  factor  is  related  to  hydrophobicity.  In  essence,  the  molecules  are 
"squeezed  out"  of  the  water  and  deposited  on  the  soil  surfaces  due  to  their  hydrophobicity 
(e.g.,  Curtis  et  al,  1986).  The  second  factor  affecting  sorption  is  the  fraction  of  solid 
organic  matter  in  the  aquifer  solids,  known  as  organic  carbon  content.  Attempts  have 
been  made  to  correlate  organic  contaminant  sorption  with  soil  organic  matter  and  a 
chemical  property  of  the  contaminant  such  as  the  octanol-water  partition  coefficient 
(Karickhoff  cr  al,  1979).  The  assumption  is  that  all  sorption  is  due  to  organic  matter 
implying  the  amount  of  sorption  is  proportional  to  the  amount  of  organic  matter. 

If  a  solute  does  not  sorb  to  the  porous  medium,  it  will  move  at  a  velocity  equal  to 
the  groimdwater  velocity  and,  by  definition,  has  a  retardation  factor  equal  to  one. 
However,  if  the  solute  does  sorb  to  the  soil,  it  will  move  at  a  velocity  equal  to  the 
groundwater  velocity  divided  by  a  retardation  factor,  R: 

R=l+|2Kd  (2-5) 

e 

where  p  =  soil  bulk  density  [Mg/L^] 

0  =  soil  porosity  [-] 

Kj  =  equilibrium  distribution  coefficient  [L^/Mj] 
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The  equilibrium  distribution  coefficient  is  a  function  of  the  soil  and  the  chemical 
type.  Numerous  attempts  have  been  made  aimed  at  correlating  with  chemical 
properties,  such  as  solubility  or  octanol-water  partitioning  coefficient  (Karickhoff,  1979; 
Chiou  et  al,  1983;  Curtis  et  al,  1986).  These  attempts  assume  that  the  sorption  process  is 
attributed  solely  to  the  soils  organic  carbon  fraction.  However,  other  studies  (e.g..  Ball 
and  Roberts,  1991a)  have  shown  that  for  soils  characterized  by  low  organic  carbon 
contents  (i.e.,  about  0.1%  or  less  w/w),  mineral  surfaces  play  a  significant  role.  In  any 
case,  the  published  correlations  are  best  used  for  a  first  approximation,  and  more  accurate 
Kd  estimates  must  be  determined  experimentally. 

2.4  Groundwater  Remediation 

Contaminated  groundwater  can  contain  an  abundance  of  organic  chemicals.  Due 
to  the  low  solubility  of  some  of  these  chemicals,  large  dilute  plumes  of  contaminant 
typically  appear  and  propagate  through  the  saturated  zone.  The  most  common  method  of 
treatment  for  a  contaminated  aquifer  is  a  pump-and-treat  system.  The  contaminated 
groundwater  is  extracted  via  a  withdrawal  well  and  treated  at  the  surface  through  a 
filtration  process  such  as  granular  activated  carbon. 

2.4.1  Effect  of  Sorption  on  Remediation 

In  an  aquifer,  one  might  expect  a  contaminant  to  move  at  the  same  speed  as  the 
groundwater;  however,  this  is  not  necessarily  the  case.  Some  contaminants  are  sorbed 
onto  the  soil  particles  thereby  reducing  the  overall  solute  flow  rate  relative  to  that  of  the 
groundwater  itself.  The  solute  is  retarded  by  a  factor  R  called  the  retardation  factor.  For 
a  successful  remediation,  R  pore  volumes  must  be  extracted  under  ideal  conditions 
(assuming  the  local  equilibrium  assumption  (LEA)  applies:  no  immobile  zones  and 
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instantaneous  equilibrium  is  achieved  between  the  water  and  solid  phases).  Thus, 
reducing  sorption  (or  R)  would  expedite  the  remediation  process. 

2.4.2  Nonequilibrium  Sorption 

Equilibrium  sorption  models  such  as  LEA  assume  that  the  mass  transfer  rate  is 
fast  relative  to  the  flow  rate  so  that  equilibrium  can  be  reached.  If  this  is  not  the  case  and 
equilibrium  is  not  attained,  a  nonequilibrium  sorption  model  is  more  appropriate. 
Unfortunately,  a  kinetic  model  causes  complications  at  two  scales.  At  the  particle  scale, 
mass  transfer  resistances  like  intraparticle  diffusion  can  markedly  decrease  the  rate  of 
contaminant  release  from  a  sorption  site  thereby  increasing  remediation  times  (e.g..  Ball 
and  Roberts,  1991a;  Harmon  and  Roberts,  1994).  At  the  layer  scale,  an  increase  in  the 
degree  of  hydraulic  conductivity  due  to  heterogeneity  (larger  scale  immobile  or  low  flow 
zones)  can  be  a  dominant  factor  in  causing  remediation  times  to  increase  (Rabideau  and 
Miller,  1994;  Kong  and  Harmon,  1996). 

2.4.3  Effect  of  Cosolvents  on  Retardation 

At  the  particle  scale,  complications  caused  by  nonequilibrium  sorption  can  be 
mitigated  by  the  use  of  cosolvents  (Brusseau  et  al,  1991).  Cosolvents  can  lower 
retardation  factors  and  increase  mass  transfer  rates.  This  thesis  will  examine  the  capacity 
of  a  cosolvent  to  lower  the  retardation  factor  of  PCE.  Before  examining  the  Research 
Objectives  and  Approach,  Cosolvency  Theory  and  Subsurface  Remediation  will  be 
reviewed  in  detail  in  the  next  chapter. 
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Chapter  5.  Cosolvency  Theory  and  Subsurface  Remediation 

The  idea  of  using  cosolvents  to  accelerate  the  remediation  of  soils  and 
groundwater  has  recently  been  considered  for  field  testing  (Augustijn  et  al,  1994). 
Although  a  relatively  new  remediation  technique  for  soil  and  groundwater,  the 
technology  was  actually  developed  in  the  petroleum  industry.  The  mechanisms  involved 
in  cosolvent  flooding  for  enhanced  oil  or  nonaqueous  phase  liquid  (NAPL)  recovery  have 
been  discussed  in  several  reviews  (Reed  et  al,  1977;  Larson  et  al,  1982;  Lake,  1983; 
Imhoff  er  al,  1995).  Other  studies  have  indicated  that  the  use  of  cosolvents  Avill  reduce 
the  retardation  factor  substantially,  thereby  vastly  reducing  the  remediation  times 
(Yalkowsky  et  al,  1981;  Rao  et  al,  1985;  &  Nkedi-Kizza  et  al,  1989). 

3.1  Methanol  as  a  Cosolvent 

Methanol  is  a  colorless,  polar  chemical  that  is  miscible  with  water.  Today, 
methanol  has  many  industrial  and  consumer  uses.  Approximately  70%  of  the  methanol 
produced  worldwide  is  used  in  chemical  syntheses  (Ullman,  1990).  The  most  readily 
apparent  use  of  methanol  is  as  an  energy  source.  The  oil  crisis  in  the  early  1970’s  set  in 
motion  the  need  to  find  alternative  sources  of  fuel.  Coupled  with  the  added  emphasis  on 
air  quality,  methanol  and  methanol-petroleum  fuel  mixtures  were  a  good  solution.  As  the 
21st  century  approaches,  the  use  of  methanol  is  increasing.  Because  of  methanol's  low 
freezing  point  and  miscibility  with  water,  it  is  sometimes  used  in  refi-igeration  systems. 
Moreover,  it  is  used  as  an  anti-ffeeze  or  an  absorption  agent  in  gas  scrubbers. 

Methanol  is  an  excellent  choice  for  a  cosolvent.  Recent  studies  (e.g.,  Imhoff  et  al, 
1995)  cite  methanol's  advantages;  (1)  it  is  relatively  inexpensive;  (2)  once  diluted  to  low 
concentrations  in  groundwater,  it  is  readily  biodegradable;  and  (3)  it  is  representative  of  a 
class  of  alcohols  that  do  not  significantly  partition  into  denser  than  water  nonaqueous 
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phase  liquids  (DNAPLS),  and  thus  will  not  enhance  downward  DNAPL  migration  by 
reducing  the  interfacial  surface  tension  of  the  DNAPL  residual  ganglia. 

3.2  Cosolvent  Theory 

The  idea  of  including  other  organic  molecules  in  water  to  affect  the  solubility  of 
organic  contaminant  in  a  solution  is  the  basis  behind  the  theory  of  cosolvency.  When  the 
organic  molecules  are  present  in  relatively  large  quantities  (s  >  10%  volume/voliune), 
they  act  as  solvent  molecules  themselves  and  partially  surround  the  solute  of  interest, 
approximately  in  proportion  to  their  volume  fraction  in  the  solution  (Yalkowsky  et  al, 
1976).  If  the  organic  molecules  are  not  present  in  these  quantities,  the  effect  on  solubility 
ranges  from  no  ejfect  (for  <  0.001  v/v)  to  a  slight  effect  (for  >  0.001  but  <  0.10). 

Yalkowsky  et  al  (1976)  postulated  that  the  excess  free  energy  of  solution  of  a 
solute  in  a  water  organic  cosolvent  mixture  should  be  a  linear  combination  of  the  solute's 
excess  free  energies  of  solution  in  each  solvent  alone.  Part  of  the  organic  solute  is 
dissolved  in  water  and  the  remainder  of  the  organic  solute  is  dissolved  in  the  organic 
cosolvent.  By  solving  for  the  excess  free  energy  of  solution  in  pure  water  and  the  excess 
free  energy  of  solution  in  the  cosolvent,  and  by  substitution,  Yalkowsky  et  al  (1976) 
showed  mathematically  that  the  log  solubility  of  a  solute  increased  linearly  as  the  fraction 
of  cosolvent  in  the  solution  mixture  increased.  This  result  was  well  supported  by 
experimental  observations  conducted  by  Yalkowsky  et  al  (1976)  using  a  series  of 
aromatic  hydrocarbons. 

3.3  Cosolvent  (Methanol)  Effects 

A  recent  study  (Imhoff  et  al,  1995)  investigated  the  following  topics:  (1)  the 
effect  of  methanol  on  various  system  parameters;  (2)  the  influence  of  methanol  on  phase 
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partitioning;  and,  (3)  the  effect  of  methanol  on  PCE  mobilization  and  rate  of  PCE 
dissolution.  Conclusions  were  made  after  conducting  a  series  of  batch-contacting  and 
generator  column  experiments. 

3.3.1  Physical  and  Chemical  Properties  of  PCE 

The  effect  of  methanol  on  various  physical  and  chemical  properties  of  PCE  were 
determined  for  a  select  number  of  methanol/water  volumetric  fractions  (0%,  20%,  40%, 
and  60%  methanol  by  volume).  The  60%  methanol  by  volume  results  demonstrated  that 
methanol  significantly  affected:  the  viscosity  of  the  aqueous  phase,  the  interfacial  tension 
between  the  aqueous  phase  and  PCE,  the  molecular  diffusion  coefficient  of  PCE,  and 
most  markedly,  the  PCE  solubility.  As  their  study  (Imhoff  et  al,  1995)  pointed  out,  an 
increased  molecular  diffusion  coefficient  and  an  increased  solubility  resulted  in  a  faster 
remediation  time. 

3.3.2  Mobilization  of  PCE 

Three  column  experiments  were  conducted  to  examine  the  effect  of  various 
methanol/water  mixtures  on  the  mobilization  of  trapped  PCE  ganglia  (Imhoff  et  al, 

1995).  Although  column  experiments  one  and  three  resulted  in  no  measurable  PCE 
mobilization,  experiment  two  demonstrated  PCE  mobility.  The  injection  of  a  60% 
methanol/water  mixture  into  the  column  increased  the  aqueous  phase  viscosity,  decreased 
the  nonaqueous-aqueous  phase  interfacial  tension,  and  increased  the  density  difference 
between  the  nonaqueous  and  aqueous  phases.  It  was  readily  apparent  in  the  analysis  of 
experiment  two  that  the  flushing  of  the  column  with  the  methanol  cosolvent  may  have 
enhanced  the  downward  migration  (mobilization)  of  DNAPL  ganglia. 
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3.3.3  Dissolution  and  Mass  Transfer  Rate  Coefficient 

Imhoff  er  al  (1995)  demonstrated  that  various  methanol/water  mixtures  had  an 
impact  on  the  rate  of  PCE  dissolution  and  the  solubility  limit  of  PCE  in  the  aqueous 
phase.  The  change  in  the  mass  transfer  rate  coefficient  with  varying  methanol  fraction 
was  due  to  the  corresponding  change  in  aqueous  phase  viscosity.  When  plotting  the 
maximum  potential  mass  flux  versus  the  methanol  fraction,  the  maximum  potential  mass 
flux  was  shown  to  increase  by  a  factor  of  30  as  the  methanol/ water  fraction  increased 
from  0%  to  60%.  The  improvement  in  the  mass  flux  was  attributed  to  the  increase  in 
aqueous  phase  PCE  solubility. 

3.3.4  Related  Research 

There  are  several  publications  on  the  effective  use  of  cosolvents  for  enhanced  in 
situ  remediation.  The  major  conclusions  from  these  are  briefly  summarized  in  this 
section. 

Nkedi-Kizza  et  al  (1985)  presented  data  in  their  study  that  clearly  demonstrated 
the  validity  of  the  solvophobic  approach  for  predicting  the  sorption  of  hydrophobic 
organic  chemicals  (HOC's)  from  binary  solvent  mixtures.  They  concluded  that  for  each 
sorbate,  the  sorption  coefficient  decreased  log  linearly  as  the  fraction  of  organic  cosolvent 
increased. 

Fu  et  al  (1986)  studied  the  hydrophobic  sorptive  behavior  of  four  aromatic  solutes 
onto  three  different  soils  with  solvent/water  mixtures.  For  a  polar  solvent,  sorption 
decreased  semi-logrithmically  with  an  increase  in  volume  fraction  of  solvent  in  the 
aqueous  phase.  This  may  have  been  the  result  of  the  solvent/water  mixture  swelling  the 
organic  carbon  associated  with  the  soil  and  thereby  increasing  solute  accessibility  to 
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organic  matter.  The  conclusion  of  the  study  was  that  the  more  hydrophobic  the  solute, 
the  greater  the  effect  of  solvent  in  solvent/water  mixtures  on  solute  solubility 
enhancement,  and  hence  the  less  the  tendency  to  sorb  onto  soil.  Thus,  the  results  of  this 
investigation  were  particularly  significant  for  those  aromatic  solutes  exhibiting  the  lowest 
aqueous  phase  solubility. 

Rao  et  al  (1990)  demonstrated  that  the  presence  of  a  completely  miscible  organic 
solvent  (CMOS)  increases  the  partially  MOS  solubility  which  in  turn,  is  reflected  in 
increased  solubility  and  decreased  sorption  of  hydrophobic  organic  chemicals  (HOC's). 
Further,  with  increased  volume  fraction  of  a  CMOS  in  a  binary  mixed  solvent,  HOC 
solubility  increased  and  sorption  decreased,  essentially  in  a  log-linear  matter.  They 
concluded  that  an  increase  in  HOC  solubility  in  the  presence  of  cosolvents  was  reflected 
by  decreased  sorption  by  soils  and  increased  mobility  of  HOC's. 

Augustijn  et  al  (1994)  demonstrated  in  laboratory  miscible  displacement 
experiments  that  the  smaller  the  retardation  factor  in  water  and  the  higher  the  cosolvent 
fraction,  the  faster  the  contaminant  was  recovered  by  solvent  flushing.  In  addition,  the 
presence  of  non-equilibrium  conditions,  soil  heterogeneity,  and  the  type  of  cosolvent 
influenced  the  time  required  to  recover  the  contaminant.  Augustijn  et  al  found  that 
solvent  flushing  was  appropriate  for  more  hydrophobic  chemicals.  Experimental  data  and 
model  simulations  showed  that  with  increasing  cosolvent  content,  the  contaminant  eluted 
at  higher  concentrations,  thus  improving  the  contaminant  recovery  efficiency.  The 
conclusion  of  the  study  was  that  high  cosolvent  fractions  were  most  effective  for  the 
elution  of  contaminants.  The  recovery  efficiency  was  expected  to  decrease  under  non¬ 
equilibrium  conditions  and  in  a  heterogeneous  soil.  The  type  of  cosolvent  and 
composition  of  the  solvent  mixture  are  design  parameters  that  can  be  used  to  optimize  the 
recovery  efficiency. 
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Harmon  et  al  (1994)  focused  on  two  aspects  of  methanol's  behavior  in  the 
subsurface  environment:  (1)  the  mixing  and  dissolution;  and,  (2)  the  mobilization  of 
sorbing  organic  contaminants.  The  studies  were  carried  out  in  batch  and  column  systems 
composed  of  a  well-characterized  sand  fraction.  All  systems  were  water-saturated.  The 
study  examined:  the  impact  of  methanol  on  the  sorption  of  the  hydrophobic  contaminant, 
benzene;  the  transport  of  methanol  and  benzene  as  solutes;  the  mixing  behavior  of 
methanol  in  a  water-saturated  soil  column;  and  the  impact  of  methanol  mobilizing 
benzene  in  a  water-saturated  soil  column.  The  results  of  the  experiments  showed  that  the 
retardation  factor  of  benzene  (i.e.,  sorption)  decreased  with  increasing  methanol  fraction. 
However,  at  higher  methanol  fractions  (>20%),  there  was  an  apparent  increase  in  the 
retardation  factor.  They  attributed  the  increase  in  retardation  to  the  density  difference 
between  water  and  methanol.  The  major  conclusions  generated  by  the  study  were:  (1) 
methanol  reduces  the  equilibrium  transport  parameter  (confirming  results  from  previous 
studies);  and  (2)  accurate  modeling  of  the  mixing  of  pure  methanol  with  water  in  porous 
media  requires  addressing  the  water-methanol  density  differences. 

3.3.5  The  Mixing  of  Miscible,  Inhomogeneous  Fluids  in  the  Subsurface 

A  major  finding  of  the  Harmon  et  al  study  (1994)  was  that  there  was  a  relatively 
poor  understanding  of  the  dynamics  of  mixing  of  miscible,  inhomogeneous  fluids  in  the 
subsurface  environment.  They  found  that  the  differences  between  methanol  and  water 
(i.e.,  density  and  viscosity)  imply  that  the  displacement  of  one  fluid  by  the  other  will  defy 
conventional  flow  and  transport  modeling  approaches,  which  treat  density  and  viscosity 
as  constants.  Currently,  research  is  under  way  by  Harmon  et  al  (1996)  to  illustrate  the 
dynamics  of  mixing  and  transport  of  miscible,  inhomogeneous  fluids  in  a  porous  medium 
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in  hopes  of  increasing  the  understanding  of,  and  design  capabilities  for,  chemically 
enhanced  remediation. 

3.4  Thesis  Objectives 

The  study  by  Harmon  et  al  (1994)  is  the  basis  for  this  thesis.  The  impact  of 
methanol  as  a  cosolvent  will  be  further  investigated;  in  this  case,  on  the  nonpolar  organic 
solute,  PCE.  Furthermore,  this  study  will  investigate  the  mixing  dynamics  of  miscible, 
inhomogeneous  fluids  (methanol  and  water)  in  the  context  of  soil  and  groundwater 
remediation  strategies. 
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Chapter  4,  Research  Objectives  and  Approach 

As  with  the  previous  study  (Harmon  et  al,  1994),  the  research  approach  for  this 
masters  thesis  focuses  on  the  study  of  a  one  dimensional,  single  aquifer  solid  type. 
However,  the  organic  contaminant  to  be  used  in  conjunction  with  methanol  is  PCE 
instead  of  benzene. 


4.1  Research  Objectives 

The  objectives  of  this  study  were  two-fold: 

(1)  To  reconfirm  the  impact  of  a  cosolvent  (i.e.,  methanol)  on  the  observed 
retardation  of  nonpolar  organic  solute  PCE  transport  in  a  water-saturated 
porous  medium. 

(2)  To  provide  evidence  that  ineffective  mixing  of  two  inhomogeneous  fluids 
(i.e.,  methanol  and  water)  can  negatively  impact  the  desired  cosolvent  effect. 

4.2  Experimental  Approach 

The  porous  media  used  in  this  study  is  Borden  sand  (-40+60  U.S.  Standard  mesh 
size,  mean  particle  diameter  =  0.33  mm).  The  experiments  performed  in  the  study  were: 

(1)  Tracer  Experiments-  Breakthrough  curves  of  trace  amoimts  of 

l4C-labeled  methanol/PCE  in  water  to  examine  their  transport  behavior. 
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(2)  Methanol  Displacement  Experiments-  A  series  of  breakthrough 
curves  for  examining  the  transport  of  various  methanol  fractions  through 
the  sand  column. 

(3)  PCE  Mobilization  Experiments-  A  series  of  breakthrough  curves  to 
observe  methanol  mobilizing  PCE  in  a  saturated  soil  column. 

4.3  Theoretical  Approach 

The  results  obtained  from  these  experiments  were  modeled  using  available 
analytical  solutions  to  the  advection-dispersion  equation.  In  the  tracer  experiments,  trace 
amounts  of  MeOH  in  water  were  run  through  the  column.  The  goal  was  to  generate 
breakthrough  curves  which  would  provide  a  basis  of  comparison  for  the  methanol 
displacement  experiments  and  PCE  mobilization  experiments. 

In  the  next  set  of  experiments,  the  goal  was  to  determine  the  advection  and 
dispersion  characteristics  of  methanol  as  a  solute  and  as  a  cosolvent.  Five  individual  tests 
were  conducted.  The  methanol  fractions  used  for  each  of  the  test  pulses  were  1%,  5%, 
10%,  20%,  and  50%.  Using  the  soil-column  lab  set-up  discussed  in  Chapter  Five,  forty 
experimental  data  points  were  measured  in  each  test  and  used  to  characterize  a  specific 
methanol  breakthrough  curve.  A  least-squares  fitting  program,  CFITM,  (van  Genuchten, 
1981)  was  used  to  find  the  Peclet  number  for  each  data  set  assuming  a  retardation  value 
of  one.  CFITM  is  a  FORTRAN  program  capable  of  solving  the  nondimensional  form  of 
the  Linear  equilibrium  adsorption  (LEA)  model  (LEA  model  discussed  in  Section  5.1). 
With  a  Peclet  number  for  each  methanol  fraction,  the  stage  would  be  set  for  the  next  set 
of  column  experiments. 
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In  the  final  set  of  experiments,  the  goal  was  to  determine  how  the  different 
methanol  fractions  would  affect  the  mobilization  of  PCE  in  the  soil  column.  Using  the 
same  experimental  set-up,  forty  experimental  data  points  were  measured  for  each 
methanol  fraction.  Using  the  Peclet  number  generated  from  the  previous  set  of 
experiments,  CFITM  was  used  to  fit  the  retardation  factor.  The  fitted  retardation  factor 
would  be  an  indication  of  methanol's  ability  to  mobilize  PCE. 
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Chapter  5,  Modeling  Approach  and  Experimental  Methods 

5.1  Modeling  Approach 

Transport  of  a  nonsorbing  solute  in  a  one-dimensional  water-saturated  porous 
medium  is  governed  by  the  advection/dispersion  equation  (e.g.,  Freeze  and  Cherry, 
1979): 


dC  —  Dhd  dC  (5"1) 

dt  dxl-  dx 

where  Djjd  =  hydrodynamic  dispersion  coefficient  [L^/T] 

Vx  =  the  average  pore  water  velocity  [L/T] 
t  =  time  [T] 

X  =  distance  [L] 

C  =  aqueous  concentration  [M/L^] 


Equation  (5-1)  assumes  steady-state  flow  and  no  interactions  between  the  aqueous  phase 
and  the  solid  phase. 

If  chemical  adsorption  is  considered,  an  additional  term  is  needed  to  account  for 
the  interaction  between  the  aqueous  phase  and  the  solid  phase  and  the  advection/ 
dispersion  equation  becomes: 

Dhd^-Vx2C-2S  (5-2) 
dt  dx^  dx  dt 

where  S  =  adsorbed  concentration  [M/M]. 


Prior  to  the  conception  of  the  CFITM  model,  significant  deviations  were  observed 
between  the  calculated  and  experimental  effluent  curves  using  equation  (5-2).  In  an 
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attempt  to  account  for  these  apparent  differences,  several  models  were  introduced  (for 
review,  see  Harmon  et  al,  1989).  The  increasing  complexity  with  which  the  immobile 
zones  were  viewed  (described  by  the  S  term)  led  to  CFITM.  (van  Genuchten,  1981). 

In  his  program  CFITM,  van  Genuchten  (1981)  considers  five  conceptual  models. 
Depending  on  the  exact  form  of  the  transport  model,  the  program  allows  up  to  five 
different  parameters  to  be  estimated  simultaneously.  The  least-squares  computer  model 
fits  these  transport  parameters  to  column  effluent  data.  The  five  transport  models  are: 
Model  A:  Linear  equilibrium  adsorption 
Model  B:  Physical  non-equilibrium 

Model  C:  Physical  non-equilibrium  in  the  presence  of  anion  exclusion 
Model  D:  Two-site  kinetic  non-equilibrium  adsorption 
Model  E:  One-site  kinetic  non-equilibrium  adsorption 
The  model  used  for  this  work  is  the  Linear  Equilibrium  Adsorption  Model  (LEA). 
With  the  LEA  model,  the  relationship  between  the  aqueous  concentration  and  the  sorbed 
concentration  is  described  by  a  linear  isotherm: 

S  =  KdC  (5-3) 

where  is  an  empirical  distribution  coefficient  [L^/MJ.  Substitution  of  equation  (5-3) 
into  equation  (5-2)  yields  the  following  transport  equation: 


R2C  =  D^-Vx2C  (5-4) 

dt  ax2  ar 


where  R  is  the  retardation  factor  [-]  equal  to: 
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(5-5) 


R=l+aKd 
0 

where  p  =  soil  bulk  density  [Mg/L^] 

0  =  soil  porosity  [-] 

K(j  ==  empirical  distribution  coefficient  [L^/Mj] 

The  program  requires  input  data  for  the  following  dimensionless  variables  in 
order  to  perform  the  analysis  of  the  effluent  data: 

(1)  Length  (Z)  =  x/L 

(2)  Concentration  (C)  =  C/Cq 

(3)  Pore  volumes  (T)  =  L 

(4)  Peclet  Number  (P)  =  v^L/  D^j) 

For  the  dimensionless  parameter  (4),  the  Peclet  number  is  defined  as  the  ratio  of 
advective  flux  to  dispersive  flux.  A  high  Peclet  number  means  advection  is  dominating 
transport  while  a  low  Peclet  number  implies  dispersion  is  dominating  transport.  Typical 
Peclet  number  values  for  natural  flow  in  groundwater  ranges  from  one  to  five;  however, 
these  values  may  be  greater  under  pumping  conditions  (Mackay  et  al,  1985). 

To  use  van  Genuchten's  model,  one  must  input  the  Peclet  number,  the  retardation 
factor,  the  length  of  the  input  (pulse),  and  the  data  points  (in  dimensionless  form).  The 
user  can  request  the  program  to  fit  one,  two  or  all  three  parameters  (Peclet  number, 
retardation  factor,  and  length  of  input).  The  fitted  parameters  are  obtained  by  means  of  a 


L  =  length  of  porous  media 


Co=  influent  concentration 
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least-squares  fit  of  the  appropriate  analytical  solution  to  column  effluent  data.  By 
iteration,  a  fitted  solution  is  obtained.  The  non-linear  least  squares  analysis  output  gives 
the  fitted  parameter  value  and  a  95%  confidence  interval  of  the  results.  Although 
relatively  old  (1981),  van  Genuchten's  program  is  not  obsolete.  The  program  is  versatile 
and  can  be  dimensioned  for  up  to  90  data  points.  The  least-squares  computer  model 
provides  a  convenient,  efficient,  and  accurate  means  of  fitting  various  transport 
parameters  to  column  effluent  data  (van  Genuchten,  1981). 

5.2  Model  Boundary  Conditions 

Analytical  solutions  of  equation  (5-4)  exist  for  several  sets  of  initial  and  boundary 
conditions.  In  van  Genuchten's  study  (1981)  for  the  LEA  model,  the  initial  condition  is: 

C(x,0)  =  Cj  (5-6) 

Two  different  conditions  can  be  applied  to  the  upper  boundary  of  the  soil  column  (x  =  0): 
a  first-type,  constant  concentration  boundary  condition  of  the  form: 

C(0,t)  =  Co  (5-7) 

or  a  second-type,  constant  flux  boundary  condition  of  the  form: 

(-D^  +  vc)  =  vCo  (5-8) 

5x  x  =  0 

where  Cq  is  the  concentration  of  the  input  solution. 

For  the  lower  boundary,  the  following  condition  is  applied: 
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5C(oo,t)  =  0 
Sx 


(5-9) 


This  condition  assumes  the  presence  of  a  semi-infinite  soil  colvimn.  An  alternative 
boundary  condition,  one  that  is  frequently  used  for  column  displacement  studies,  is  that 
of  a  zero  concentration  gradient  at  the  lower  end  of  the  column: 

2C(L,t)  =  0  (5-10) 

dx 

where  L  is  the  column  length. 

In  van  Genuchten's  study  (1981),  equations  (5-9)  and  (5-10)  are  both  used. 
However,  only  the  solutions  for  a  semi-infinite  medium  are  included  in  the  curve-fitting 
program.  Because  of  the  relatively  small  influence  of  the  imposed  mathematical 
boundary  conditions,  the  solutions  for  a  semi-infinite  system  provide  close 
approximations  for  those  applicable  to  a  physically  well-defined  finite  system,  especially 
for  not  too  short  laboratory  soil  columns. 

5.5  Experimental  Methods 

Each  of  the  experiments  summarized  in  Chapter  4.3  is  described  in  detail  in  the 
following  section. 

5.3.1  One-Dimensional  Column  Set-up 

The  set-up  of  the  soil  column  experiments  is  shown  schematically  in  Figure  5-1. 
The  controlled  flow  of  methanol  and  water  was  driven  through  the  soil  column  using  a 
quartemary  high  performance  liquid  chromatography  (HPLC)  pump  (Paramus,  NJ)  in 
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SYRINGE  PUMP 


Figure  5-1.  Schematic  diagram  of  column  study  experimental  set-up 
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order  to  deliver  precise  mixtures  of  methanol  and  water.  In  this  study,  a  25-cm  steel 
column  (2.06  cm  inner  diameter)  was  used,  packed  with  -40+60  (U.S.  Std.  Mesh  Size) 
Borden  sand,  which  had  been  examined  in  particle  characterization  and  organic 
contaminant  sorption  and  desorption  studies  with  PCE  (Ball  et  al,  1990;  Ball  and 
Roberts,  1991a;  1991b;  Harmon  and  Roberts,  1994).  The  Borden  sand  has:  a  mean  grain 
diameter  of  0.33  mm;  a  density  of  2.71  g/cm^;  and  an  organic  carbon  content  of 
0.00023%  (Ball  et  al,  1990).  The  porosity  was  approximately  0.35.  The  l^C-labeled 
solute,  either  methanol  or  PCE,  was  injected  via  a  gastight  syringe  mounted  on  a  syringe 
pump  (Sage  Model  34 IB,  Boston,  MA).  A  three-way,  electronically  actuated  valve 
(Valeo,  Houston,  TX),  which  allowed  precise  switching  between  the  HPLC  pump  and  the 
syringe  pump,  was  used  for  inputting  the  I'lC-labeled  solute  pulse.  The  data  points  for 
the  effluent  breakthrough  curves  were  captured  in  a  series  of  vials  containing  liquid 
scintillation  cocktail  fluid.  All  fittings  were  composed  of  glass.  Teflon,  and  stainless 
steel  to  avoid  interactions  between  the  organic  contaminants  and  the  apparatus. 

5.3.2  Analytical  Methods 

5.3.2.1  Liquid  Scintillation  Counting 

Liquid  scintillation  vials  typically  contained  5  mL  of  liquid  scintillation  cocktail 
(Insta-gel,  Packard  Instrument  Co.,  Downer's  Grove,  IL).  The  samples  were  counted  with 
a  Packard  Tricarb  Model  4530  scintillation  counter  (Packard  Instrument  Co.).  All 
samples  were  counted  two  times  for  20  minutes  each.  Measured  counts  per  minute  were 
converted  to  actual  disintegrations  per  minute  by  the  external  standard  method,  with 
quenching  efficiencies  estimated  by  a  channels  ratio  method  incorporated  into  the  system 
software.  Quenching  was  less  than  5%  for  all  samples  and  in  a  spectral  region  where 
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quenching  was  relatively  insensitive  to  changes  in  solution  matrix.  Each  set  of  vials  to  be 
counted  was  accompanied  by  a  blank  sample  (containing  5  mL  of  scintillation  cocktail), 
used  to  subtract  out  background  activity. 

5.3.2.2  Gas  Chromatography 

Prior  to  each  spiking  session,  the  overall  (labeled  and  non-labeled)  concentration 
of  the  spiking  solution  was  determined  by  comparison  with  standard  solutions  using  gas 
chromatography.  Samples  of  the  spiking  solution  were  prepared  by  extracting 
approximately  3  mL  of  the  spiking  solution  with  0.5  mL  pentane  in  2  mL  vials  equipped 
with  teflon-lined  silicon  septa  and  plastic  screw-top  caps  (Alltech  Assoc.,  Deerfield,  IL). 
Samples  and  standards  were  also  spiked  with  an  internal  standard  (1-chloro,  2- 
bromopropane).  The  PCE  and  internal  standard  were  extracted  into  the  pentane  during  30 
minutes  of  agitation  on  a  shaker  table. 

The  gas  chromatograph,  a  Hewlett-Packard  5 890 A,  was  equipped  with  a  capillary 
column.  The  column  was  followed  by  a  ^^Ni  electron  capture  detector  (ECD).  The 
splitting  valve  was  used  on  occasion  to  test  the  purity  of  new  purchases  of 
radiochemicals.  The  valve  provided  the  capability  of  diverting  most  of  the  flow  (99%) 
away  from  the  detector  in  these  instances,  when  a  relatively  concentrated  and  radioactive 
injection  was  applied.  For  quantitative  purposes,  2-7  pL  injections  of  the  pentane  were 
applied,  and  much  less  of  the  flow  (<5%)  was  diverted  from  the  detector.  The  carrier  and 
makeup  gases  were  helium  and  Argon/Methane,  respectively.  The  system  operated  at 
150  degrees  Celsius,  with  the  detector  slightly  higher-  160  degrees  Celsius. 
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5.3.2.3  Radiochemical  Purity 

Prior  to  spiking  procedures,  spiking  solutions  were  routinely  analyzed  for  volatile 
impurities  via  the  gas  chromatography  procedures  described  in  the  preceding  section.  No 
extraneous  peaks  were  detected  in  analyses  with  the  PCE  spiking  solutions. 

Spiking  solutions  were  routinely  tested  for  I'^C-nonvolatile  and  ^^C-C02  fractions 
using  an  acidd?ase  purging  procedure  (Lanzarone  and  McCarty,  1988).  The  nonvolatile 
impurities,  which  were  not  identified,  were  found  to  represent  approximately  1-2%  of  the 
activity,  for  the  PCE.  This  fraction  appeared  to  remain  in  the  aqueous  phase  (i.e.,  did  not 
sorb)  in  samples  for  which  the  A-B-N  test  was  performed  prior  to  purging.  No  significant 
fraction  of  14C-C02  was  detected  in  the  PCE  spiking  solutions.  The  equilibrium 
measurement  of  the  sample  supernatant  was  adjusted  downward  to  account  for  the 
nonvolatile  fraction. 

5.3.3  Column  Experiment  Protocol 

5.3.3.1  Preliminary  Preparation 

Calibrate  the  scintillation  counting  machine  (Beckman  LS600,  Fullerton,  CA).  At 
the  quartemary  HPLC  pump,  fill  the  water  tank  with  distilled  water  (reservoir  B)  and  fill 
the  methanol  tank  with  methanol  (reservoir  C).  Degas  the  water  and  methanol  tanks  with 
helium  for  a  1/2  hour.  Turn  on  helium  gas  at  main  tank  valve  and  the  brass  toggle  valve 
that  branches  to  the  Teflon  tubing  to  the  quartemary  HPLC  pump.  Turn  on  the  helium 
flow  at  the  pump  and  get  a  flow  of  helium  bubbles  in  the  reservoirs.  During  degassing, 
continue  lab  preparation.  Label  vials  (1  through  40)  and  fill  with  5  mL  of  scintillation 
cocktail.  Record  the  weight  of  each  of  the  40  vials  before  the  experiment  begins. 
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5.3.3.2  Solute  Preparation 

Retrieve  the  l^C-labeled  PCE  or  the  l^C-labeled  methanol  sample  from  the 
refrigerator  and  allow  the  sample  to  equilibrate  at  ambient  temperature.  A  glass  syringe 
(Hamilton,  Reno,  NV)  was  used  to  load  10  pL  of  PCE  (or  50  pL  of  methanol)  into  the  10 
mL  syringe  pump  (filled  with  degassed  water  or  precise  methanol/water  mixture). 
Preliminary  tests  indicated  that  these  amounts  were  sufficient  solute  inputs.  All  syringes 
should  be  cleaned  using  a  double-acetone  wash  procedure  after  each  use.  Mix  the  1  Re¬ 
labeled  solute  input  (10  mL  syringe  pump  spiked  with  radio-labeled  PCE  or  radio- 
labeled  methanol)  over  a  magnetic  mixer.  Replace  sample  caps  with  new  teflon  caps. 

Find  Cq  (initial  concentration  in  syringe  pump).  Using  three  vials  labeled  A,  B, 

and  C;  manually,  via  a  syringe,  load  5  pL  of  the  radio-labeled  10  mL  sample  into  each  of 
the  three  vials.  Cap  and  shake  immediately.  This  procedure  is  performed  three  times  (A, 
B,  &  C)  to  account  for  possible  non-uniform  mixing  in  the  syringe  pump.  Clean  the 
syringe  with  the  double-acetone  wash  procedure  after  each  sample:  A,  B,  &  C.  Using 
the  counting  machine,  obtain  a  "quick  dpm  (disintegrations  per  minute)  count"  on  A,  B, 

&  C  vials  to  verify  a  uniform  C^.  A  "quick  dpm  count"  means  the  counting  machine, 

using  a  smaller  sample  time,  will  give  a  dpm  count  that  still  maintains  a  relative  error  that 
is  less  than  2%.  The  initial  concentration  is  calculated  by  finding  the  average  of  the 
concentrations  of  A,  B,  and  C. 

When  the  three  counts  are  within  5%  of  each  other,  proceed  with  the  experiment. 
Set  the  quartemary  HPLC  pump  to  the  desired  flow  rate  (0.745  ml/min)  and  desired 
methanol/water  mixture  (0%,  1%,  5%,  10%,  20%,  or  50%).  Load  empty  vials  1  through 
40  in  the  auto  sampler  (Isis  Autosampler,  Lincoln,  NE).  Load  about  three  rows  of  empty 
test  vials  in  the  auto  sampler.  Set  sample  and  wash  times  on  the  auto  sampler  (i.e.,  20 
seconds  sample  interval  and  one  minute  35  seconds  wash  time). 
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5.3.3.3  Lab 

Record  start  time  of  the  experiment.  Actuate  the  syringe  pump  to  load  the  step 
input  into  the  column.  When  the  syringe  is  completely  empty,  turn-off  syringe  pump, 
switch  the  actuator,  and  engage  the  quartemary  HPLC  pump.  Note  the  time  the  pulse 
was  stopped.  This  is  the  length  of  the  step  input.  Note  the  time  when  the  supply  of  test 
vials  is  exhausted  and  the  time  when  the  first  recorded  sample  is  taken. 

While  the  experiment  is  in  progress,  monitor  the  quartemary  HPLC  pump  for 
stoppages.  Note  the  time  of  any  mishaps  during  the  experiment  run.  Continue  to  sample 
the  remaining  39  vials.  Cap  and  shake  each  vial  after  every  sample  is  taken.  Rinse 
sample  tube  after  each  sample  vial  is  taken  using  distilled  water.  At  the  end  of  the 
experiment,  note  the  time  of  completion.  After  the  last  sample,  turn  off  the  auto  sampler 
while  the  sampling  rod  is  in  the  washing  tube. 

53.3.4  Post-Lab 

Weigh  each  of  the  40  sample  vials.  Load  the  A,  B,  C  vials  and  the  40  sample 
vials  into  the  scintillation  counter.  Program  the  counting  machine  to  count  the  dpm's  (ten 
minute  counting  time,  twice  each  vial).  The  packed  column  must  be  flushed  completely 
with  degassed  water  before  another  experiment  is  run  (preferably  the  next  day).  When 
flushing  the  column,  ensure  the  packed  column  is  in  the  vertical  position.  Turn  off 
helium  sources  at  the  HPLC  pump  and  the  helium  tank.  Turn  off  all  machines.  Clean-up 
covmters  and  police  the  area. 

5.4  Data  Reduction 

The  data  from  the  scintillation  counter  will  give  the  "dpm  coimt"  for  each  of  the 
sample  vials.  The  "dpm  count"  is  a  mass  measurement.  To  be  used  in  the  CFITM  model. 


33 


the  mass  measurement  must  be  converted  to  a  concentration  measurement.  Each  "dpm 
count"  will  be  divided  by  the  difference  between  the  final  and  initial  weights  of  each 
corresponding  sample  vial.  The  difference  between  the  final  and  initial  weight  of  each 
sample  vial  is  in  the  units  of  grams.  To  convert  this  difference  measurement  to  mLs, 
divide  by  the  density  of  water  (=  1  g/mL).  The  "dpm/mL"  measurement  is  now  in 
concentration  (aqueous)  units.  All  data  reduction  spreadsheets  are  included  in  Appendix 
A,  Column  Run  Spreadsheets. 

5.5  CFITM  Model  Input 

The  data  input  for  the  model  for  each  sample  vial  is  equivalent  "pore  volume"  [-], 
and  corresponding  sample  concentration  (in  [dpm/mL];  as  calculated  in  section  5.4). 
CFITM  provides  a  least-squares  fit  for  each  experimental  data  point.  This  output,  when 
connected,  is  the  fitted  curve  for  the  column  effluent  data. 

The  CFITM  data  input  files  for  each  of  the  column  runs  is  shown  in  Appendix  B, 
CFITM  Data  Input. 

5.6  Expected  Shape  of  Breakthrough  Curves 

The  processes  affecting  solute  transport  in  the  saturated  zone  are  advection  and 
hydrodynamic  dispersion.  Advection  would  carry  the  dissolved  chemical  with  the 
groundwater  flow.  In  an  ideal  porous  media,  it  would  produce  close  to  "plug  flow" 
behavior.  Hydrodynamic  dispersion  would  cause  the  dissolved  chemical  to  spread 
throughout  the  porous  media.  The  curves  in  Figure  5-2  show  the  shape  of  the  expected 
breakthrough  curve  given  a  square  pulse  input  of  solute. 
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square  input  pulse  output  breakthrough  curve 


Figure  5-2.  Expected  output  breakthrough  curve  given  a  square  input  pulse. 

In  Figure  5-2,  tg  is  the  corresponding  time  for  the  center  mass  of  the  square  input  pulse. 
The  degree  of  "spreading"  in  the  output  breakthrough  curve  varies  according  to  the 
amount  of  hydrodynamic  dispersion.  As  shown  on  the  time-scale  of  the  output 
breakthrough  curve,  the  arrival  time  of  the  peak  is  a  function  of  the  length  of  and  the  rate 
at  which  the  input  pulse  travels. 


Chapter  6.  Results  and  Discussion 


6.1  Overview 

Specific  results  are  described  in  the  following  sections.  Section  6.2  presents 
column  breakthrough  data  for  the  miscible  displacement  tests;  that  is,  employing 
methanol  or  PCE  as  single  solutes.  Section  6.3  presents  the  results  for  the  methanol  slug 
tests  in  which  cosolvent  levels  of  methanol  were  used  to  displace  water  from  the  column. 
Section  6.4  presents  the  PCE  breakthrough  curves  subject  to  MeOH  cosolvent  levels. 
Throughout  Chapter  6,  the  breakthrough  data  are  plotted  along  with  CFITM  model  output 
derived  from  the  curve  fitting  protocols  discussed  in  Section  5.5.  Output  from  the 
CFITM  program  for  each  column  run  is  included  in  Appendix  C,  CFITM  Data  Output. 

The  breakthrough  curves  for  each  of  the  column  runs  are  shown  in  the  following 
sections.  As  discussed  previously  in  Section  5.6,  all  of  the  breakthrough  curves  are 
similar  in  shape  to  the  output  curves  shown  in  Figure  5-2. 

6.2  Miscible  Displacement  Results:  Methanol  and  PCE 

The  first  two  sets  of  column  experiments  were  tracer  studies  employing  dilute 
amounts  of  either  l^C-labeled  MeOH  or  l^C-labeled  PCE.  In  each  case,  the  input  pulse 
contained  no  methanol  (0%  methanol  fraction).  The  purpose  of  the  tracers  was  to  provide 
a  basis  for  comparison  of  the  methanol  displacement  experiments  and  the  PCE 
mobilization  experiments.  Figure  6-1  shows  the  tracer  run  for  the  l^C-labeled  MeOH. 
Figure  6-2  shows  the  tracer  run  for  the  I'^C-labeled  PCE. 
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Methanol  Tracer  Run 


0.8  T 


Figure  6-1 .  Methanol  tracer  breakthrough  curve  for  a  0.745  mL/min  flowrate. 


PCE  Tracer  Run 


Figure  6-2.  PCE  tracer  breakthrough  curve  for  a  0.745  mL/min  flowrate. 
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The  solid  squares  on  Figure  6-1  and  Figure  6-2  represent  the  experimental  data.  The 
solid  line  represents  CFITM's  best-fitting  curve,  using  a  Peclet  number  equal  to  159.  A 
summary  of  the  methanol  and  PCE  tracer  tests  are  shovvn  below  in  Table  6-1 . 


Column 

Methanol 

Flow 

Peclet 

Retardation 

95%  Confidence  Limit 

Run 

(%) 

(mL/min) 

Number 

Factor 

Lower 

Upper 

Tracer  MeOH 

0 

0.745 

159 

1.00 

112 

206 

Tracer  PCE 

0 

0.745 

159 

1.20 

1.18 

1.22 

Table  6-1 .  Summary  of  Methanol  and  PCE  Tracer  Runs. 


Note-  Numbers  in  bold  italics  indicate  fitted  data. 

For  each  of  the  methanol  tracer  runs,  the  breakthrough  curve  begins  to  appear  at 
approximately  one  pore  volume.  The  expectation  is  that  the  center  of  the  breakthrough 
curve  should  elute  at  one  pore  volume  plus  half  of  the  total  length  of  the  input  pulse  (also 
in  pore  volumes).  The  center  of  the  methanol  tracer  breakthrough  curve  occurs,  for  each 
tracer  run,  at  approximately  1.5  pore  volumes.  The  peaks  arrive  slightly  later  than 
expected;  however,  considering  experimental  error,  the  tracer  curves  show  that  the 
methanol  was  transported  with  the  pore  water,  and  not  retarded. 

The  experimental  tracer  runs  exhibited  predictable  behavior.  The  0.745  mL/min 
flow  rate  used  in  each  of  the  colunm  runs  was  a  relatively  fast  flow  rate  as  shown  by  the 
shape  of  the  breakthrough  curves.  The  sharpness  of  the  front  and  elution  portions  of  the 
curves  were  indicative  of  advection  dominated  transport  characteristic  of  a  high  Peclet 
number. 
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6.3  MeOH  Displacement  Experiments 

The  objective  of  the  second  set  of  column  experiments  was  to  characterize  the 
advection  and  dispersion  behavior  of  methanol  as  a  cosolvent.  More  specifically,  the 
purpose  was  to  characterize  the  extent  that  methanol  displacing  water  behaved 
immiscibly,  due  to  density  and  viscosity  differences.  The  methanol  fi'actions  tested  (1%, 
5%,  10%,  20%,  and  50%  v/v)  exhibited  similar  behavior.  Figures  6-3  through  6-7  show 
the  methanol  breakthrough  curves  for  the  1%,  5%,  10%,  20%,  and  50%  methanol 
fractions  respectively. 


MeOH  Slug  (1%  MeOH/99%  Water) 


0.8  T 


Figure  6-3.  1%  MeOH  slug  breakthrough  at  a  0.745  mL/min  flowrate. 
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MeOH  Slug  (5%  MeOH/95%  Water) 


Figure  6-4.  5%  MeOH  slug  breakthrough  at  a  0.745  mL/min  flowrate. 


MeOH  Slug  (10%  MeOH/90%  Water) 


Figure  6-5.  10%  MeOH  slug  breakthrough  at  a  0.745  mL/min  flowrate. 
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MeOH  Slug  (20%  MeOH/80%  Water) 


Figure  6-6.  20%  MeOH  slug  breakthrough  at  a  0.745  mL/min  flowrate. 


MeOH  Slug  (50%  MeOH/50%  Water) 


Figure  6-7.  50%  MeOH  slug  breakthrough  at  a  0.745  mL/min  flowrate. 


The  simulation  providing  the  best-fit  for  the  data  was  characterized  by  a  Peclet  number  of 
328,  the  resulting  Peclet  number  for  the  20%  MeOH  fraction.  Each  of  the  column  runs 
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had  similar  square  mass  inputs  (pulses).  When  calculating  the  area  below  each  of  the 
output  curves,  it  was  estimated  that  at  least  93%  of  the  mass  input  was  conserved  for  each 
colvunn  run  (Appendix  A,  Column  Run  Spreadsheets). 

Each  of  the  column  runs  produced  breakthrough  curves  that  were  characteristic  of 
advection  dominated  transport.  Analysis  of  the  breakthrough  curves  shows  that  there  was 
no  evidence  of  unstable  fronts  nor  was  there  evidence  of  "over-riding"  or  "wedging"  of 
methanol  over  water.  Table  6-2  summarizes  the  results  of  the  MeOH  slug  experiments. 


Column 

Run 

Methanol 

(%) 

Flow 

(mL/min) 

Peclet 

Number 

Retardation 

Factor 

95%  Confidence  Limit 

Lower 

Upper 

MeOH  Slug 

1 

0.745 

153 

1.00 

106 

200 

MeOH  Slug 

5 

0.745 

305 

1.00 

224 

385 

MeOH  Slug 

10 

0.745 

276 

1.00 

199 

352 

MeOH  Slug 

20 

0.745 

328 

1.00 

293 

364 

MeOH  Slug 

50 

0.745 

221 

1.00 

190 

253 

Table  6-2.  Summary  of  MeOH  Slug  Experiments. 


Note-  Numbers  in  bold  italics  indicate  fitted  data. 

The  MeOH  slug  run  at  1%  methanol  by  volume  was  characterized  by  a  Peclet  number 
similar  to  that  of  the  MeOH  tracer  run.  However,  the  methanol  runs  with  greater 
methanol  percentages  by  volume  were  characterized  by  significantly  larger  Peclet 
numbers.  The  sharpness  of  the  corresponding  curve  fronts  were  also  indicative  of  the 
larger  Peclet  numbers.  Each  of  the  methanol  runs  exhibited  early  breakthrough,  however. 
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this  was  due  to  improperly  flushing  the  radio-labeled  MeOH  from  the  previous  column 
experiments.  An  adjusted  C/Cq  value  was  calculated  to  correct  for  this  inaccuracy. 

6.4  PCE  Mobilization  Breakthrough  Curves 

The  final  set  of  experiments  determined  the  effect  of  different  methanol  fractions 
on  the  mobilization  of  PCE  in  the  soil  column.  Again,  the  methanol  fractions  tested  (1%, 
5%,  10%,  20%,  and  50%  v/v)  exhibited  similar  behavior.  From  the  output  curves,  one 
would  expect  that  with  increasing  methanol  fraction,  the  curves  would  exhibit  earlier 
breakthrough,  indicative  of  a  decrease  of  the  retardation  factor  of  PCE.  Figures  6-8 
through  6-12  show  the  PCE  mobilization  breakthrough  curves  for  each  of  the  methanol 
fractions. 


PCE  Mobil.  (1%  MeOH/99%  Water) 


Figure  6-8.  PCE  mobilization  breakthrough  (1%  MeOH)  at  a  0.745  mL/min  flowrate. 
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PCE  Mobil.  (5%  MeOH/95%  Water) 


Figure  6-9.  PCE  mobilization  breakthrough  (5%  MeOH)  at  a  0.745  mL/min  flowrate. 


PCE  Mobil.  (10%  MeOH/90%  Water) 


Figure  6-10.  PCE  mobilization  breakthrough  (10%  MeOH)  at  a  0.745  mL/min  flowrate. 
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PCE  Mobil.  {20%  MeOH/80%  Water) 


Figure  6-11.  PCE  mobilization  breakthrough  (20%  MeOH)  at  a  0.745  mL/min  flowrate. 


PCE  Mobil.  (50%  MeOH/50%  Water) 


Figure  6-12.  PCE  mobilization  breakthrough  (50%  MeOH)  at  a  0.745  mL/min  flowrate. 


All  of  the  PCE  retardation  results  (1%,  5%,  10%,  20%,  and  50%)  were  similar.  As  with 
the  MeOH  slug  runs,  each  of  the  PCE  mobilization  column  runs  produced  breakthrough 
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curves  that  were  characteristic  of  advection  dominated  transport.  The  simulation 
providing  the  best-fit  for  the  data  was  characterized  by  a  retardation  factor  of  1 .22, 
corresponding  to  the  20%  MeOH  fraction.  Each  of  the  column  runs  had  similar  square 
mass  inputs  (pulses).  When  calculating  the  area  below  each  of  the  output  curves,  it  was 
estimated  that  at  least  97%  of  the  mass  input  was  conserved  for  each  column  run.  Table 
6-3  summarizes  the  PCE  retardation  results. 


Column 

Run 

Methanol 

(%) 

Flow 

(mL/min) 

Peclet 

Number 

Retardation 

Factor 

95%  Confidence  Limit 

Lower 

Upper 

PCE  Input 

1 

0.745 

153 

1.22 

1.21 

1.24 

PCE  Input 

5 

0.745 

305 

1.25 

1.25 

1.26 

PCE  Input 

10 

0.745 

276 

1.25 

1.24 

1.26 

PCE  Input 

20 

0.745 

328 

1.23 

1.22 

1.24 

PCE  Input 

50 

0.745 

221 

1.24 

1.23 

1.25 

Table  6-3.  Summary  of  PCE  Retardation  Experiments. 


Note-  Numbers  in  bold  italics  indicate  fitted  data. 

In  general,  the  PCE  breakthrough  curves  were  symmetrical  and  weakly  retarded.  The 
PCE  runs  with  20%  methanol  and  50%  methanol  each  exhibited  slightly  early 
breakthrough  and  slight  tailing. 

Van  Genuchten's  Least  Squares  Program  determined  the  best-fitting  retardation 
factor  for  each  of  the  observed  breakthrough  curves.  All  of  the  resulting  retardation 
factors  are  plotted  in  Figure  6-13  as  a  function  of  the  methanol  fraction. 
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Figure  6-13.  Retardation  Factors  with  respect  to  methanol  fractions. 


The  results  do  not  indicate  any  trends.  The  retardation  factor  only  varied  between  1 .20 
and  1.25.  Essentially,  there  was  no  significant  difference  in  retardation  factor  with 
increasing  methanol  fraction. 

The  PCE  mobilization  experiments  failed  to  demonstrate  the  hypothesized  results. 
The  degree  of  retardation  of  PCE  did  not  decrease  log-linearly  with  increasing  methanol 
fraction.  However,  the  results  of  this  study  confirmed  the  results  of  previous  studies 
(Imhoff,  1995).  The  MeOH  slug  experiments  showed  that  with  an  increasing  methanol 
fraction,  there  was  an  increase  in  the  Peclet  number. 
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Chapter  7.  Summary  and  Conclusions 

The  goal  of  this  research  was  to  test  the  ability  of  methanol  as  an  effective 
cosolvent.  This  study  investigated  the  transport  and  dissolution  of  methanol  in  a  water- 
saturated  porous  medium  and  the  effects  of  methanol  on  the  mobilization  of  PCE.  1  -D 
column  experiments  were  used  to  examine;  the  transport  behavior  of  trace  amounts  of 
14C-labeled  methanol/PCE  in  water;  the  transport  of  various  methanol  fractions  through  a 
sand  column;  and,  methanol  mobilizing  PCE  in  a  saturated  soil  column. 

The  one-dimensional  column  apparatus  devised  for  this  work  performed 
adequately.  An  average  of  95%  of  the  mass  input  through  the  column  was  calculated  in 
the  mass  output  indicating  that  a  majority  of  the  mass  was  conserved.  The  tracer  and  slug 
breakthrough  curves  of  the  methanol  were  described  well  by  the  advection/dispersion 
equation,  and  showed  no  indications  of  density  and  viscosity  difference  effects. 

For  the  1-D  laboratory  column  at  the  flowrate  tested  (0.745  mL/min),  methanol 
failed  to  produce  the  hypothesized  cosolvent  effect  in  mobilizing  PCE.  Two  possible 
reasons  for  this  failure  exist:  (1)  The  PCE  sorbed  too  weakly  to  the  porous  medium  and 
any  reductions  in  the  sorption  were  difficult  to  observe  given  the  experimental 
uncertainty  produced  by  the  methods;  (2)  The  MeOH  front  by-passed  much  of  the  porous 
medium,  thereby  failing  to  mix  sufficiently  in  a  timely  manner.  Given  the  apparent 
stability  of  the  methanol  slug  breakthrough  curves,  the  first  reason  appears  to  be  the  more 
likely  of  the  two. 

Further  studies  on  flow  conditions  in  a  contaminated  region  need  to  be 
investigated.  Perhaps  the  0.745  mL/min  flow  rate  used  in  this  study  was  too  fast  to 
adequately  demonstrate  the  impact  of  methanol  on  sorption  rates.  Although  the  outcome 
of  this  research  was  not  completely  as  expected,  the  results  do  provide  motivation  for 
future  work  to  address  these  flow  and  transport  issues. 
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Appendix  A.  Column  Run  Spreadsheets 
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Explanation  of  Terms  for  Spreadsheets- 


Type  of  column  run.  First  number  refers  to  MeOH  fraction;  second  number 
refers  to  percent  water. 

flow:  flow  rate  at  which  the  experiment  was  conducted. 

time:  the  time  at  which  the  sample  was  completed.  For  each  of  the  column  runs, 

the  sample  time  in  the  vial  was  20  seconds,  the  wash  time  in  the  waste  stream  was 
100  seconds,  and  the  time  in  between  sampling  and  washing  was  3  seconds.  Total 
time  between  sampling  vials  was  123  seconds  (2  minutes,  3  seconds). 

conv.:  time  conversion  to  decimals. 

pore  vol.:  time  conversion  to  equivalent,  unitless  "pore  volumes." 

dpm:  "disintegrations  per  minute"  as  recorded  by  the  scintillation  counter,  (dpm  = 
concentration  [C]). 

diff:  Difference  in  weight  [g]  between  the  sample  vial  initially  and  the  sample  vial  after 
sampling. 

dpmydiff:  Self-explanatory.  Assuming  the  density  of  water  is  1  g  per  cm^,  [dpm/mL]. 

Co:  Initial  concentration  of  input  pulse  [dpm/mL]. 

C/Co:  Dimensionless  parameter  [sample  concentration/initial  concentration]. 

Adj  C/Co:  Adjusted  C/Co  value  to  accoimt  for  inaccuracies  caused  by  insufficient 
column  flushing  [-]. 

M:  Mass  of  radio-labeled  MeOH  or  PCE  in  the  sample  vial  [dpm] . 

Mo:  Initial  mass  of  input  pulse  [dpm]. 

M/Mo:  Dimensionless  parameter  [sample  mass/total  initial  mass]. 

CFIT:  Fitted  value  by  van  Genuchten's  model  based  on  experimental  output  [dpm/mL]. 
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Appendix  B,  CFITM  Input  Files 
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1 

2  1  25  40 

Methanol  Run:  0%  Methanol/100%  Water 
PECLET  RF  PULSE 

5.0  1.0  0.319 

110 
0.400  0.0000 

0.451  0.0000 

0.503  0.0000 

0.554  0.0000 

0.605  0.0000 

0.656  0.0000 

0.708  0.0000 

0.759  0.0000 

0.810  0.0000 
0.861  0.0000 
0.913  0.0000 

0.964  0.0028 

1.015  0.0111 

1.066  0.0184 

1.118  0.0320 

1.169  0.0563 

1.220  0.0911 

1.271  0.1527 

1.323  0.2566 

1.374  0.4205 

1.425  0.5675 

1.476  0.6260 

1.528  0.6289 

1.579  0.5736 

1.630  0.4732 

1.681  0.3275 

1.733  0.1455 

1.784  0.0284 

1.835  0.0000 

1.886  0.0000 

1.938  0.0000 

1.989  0.0000 

2.040  0.0000 

2.091  0.0000 

2.143  0.0000 

2.194  0.0000 

2.245  0.0000 

2.296  0.0000 

2.348  0.0000 

2.399  0.0000 

C$stop 
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1 

2  .1  25  40 

Methanol  Run:  1%  Methanol/99%  Water 
PECLET  RF  PULSE 

5.0  1.0  0 . 326 

110 
0.431  0.0000 

0.483  0.0000 

0.534  0.0000 

0.585  0.0000 

0.636  0.0000 

0.688  0.0000 

0.739  0.0000 

0.790  0.0000 

0.841  0.0000 

0.893  0.0000 

0.944  0.0000 

0.995  0.0046 

1.046  0.0136 

1.098  0.0245 

1.149  0.0436 

1.200  0.0820 
1.251  0.1352 

1.303  0.2179 

1.354  0.3383 

1.405  0.4895 

1.456  0.6194 

1.508  0.6339 

1.559  0.6113 

1.610  0.5364 

1.661  0.4029 

1.713  0.2162 

1.764  0.0678 

1.815  0.0000 

1.866  0.0000 

1.918  0.0000 

1.969  0.0000 

2.020  0.0000 
2.071  0.0000 

2.123  -  0.0000 
2.174  0.0000 

2.225  0.0000 

2.276  0.0000 

2.328  0.0000 

2.379  0.0000 

2.430  0.0000 

C$stop 
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1 

2  1  25  40 

Methanol  Run:  5%  Methanol/95%  Water 
PECLET  RF  PULSE 

5.0  1.0  0 . 326 

110 
0.456  0.0000 

0.508  0.0000 

0.559  0.0000 

0.610  0.0000 
0.661  0.0000 
0.713  0.0000 

0.764  0.0000 

0.815  0.0000 

0.866  0.0000 
0.918  0.0000 

0.969  0.0000 

1.020  0.0000 

1.071  0.0024 

1.123  0.0108 

1.174  0.0318 

1.225  0.0926 

1.276  0.2447 

1.328  0.4720 

1.379  0.6573 

1.430  0.7483 

1.481  0.7455 

1.533  0.7510 

1.584  0.6581 

1.635  0.4568 

1.686  0.1993 

1.738  0.0302 

1.789  0.0000 

1.840  0.0000 

1.891  0.0000 

1.943  0.0000 

1.994  0.0000 

2.045  0.0000 

2.096  0.0000 

2 . 148  ■  0 . 0000 
2.199  0.0000 

2.250  0.0000 

2.301  0.0000 

2.353  0.0000 

2.404  0.0000 

2.455  0.0000 

C$stop 


66 


1 

2  .  1 

25  40 

Methanol 

Run:  10% 

PECLET 

RF 

5.0 

1.0 

1  1 

0 

0.394 

0.0000 

0.445 

0.0000 

0.496 

0.0000 

0.548 

0.0000 

0.599 

0.0000 

0.650 

0.0000 

0.701 

0.0000 

0.753 

0.0000 

0 . 804 

0.0000 

0.855 

0.0000 

0.906 

0.0000 

0.958 

0.0000 

1.009 

0.0000 

1.060 

0.0018 

1.111 

0.0041 

1.163 

0.0353 

1.214 

0.1728 

1.265 

0.4012 

1.316 

0.6129 

1.368 

0.7223 

1.419 

0.8058 

1.470 

0.7435 

1.521 

0,6161 

1.573 

0.5071 

1.624 

0.3088 

1.675 

0.0880 

1.726 

0.0000 

1.778 

0.0000 

1.829 

0.0000 

1.880 

0.0000 

1.931 

0.0000 

1.983 

0.0000 

2.034 

0.0000 

2.085 

-  0.0000 

2 . 136 

0.0000 

2.188 

0.0000 

2.239 

0.0000 

2.290 

0.0000 

2.341 

0.0000 

2.393 

0.0000 

C$stop 


Methanol/90%  Water 
PULSE 
0.328 


67 


1 

2  .1  25  40 

Methanol  Run:  20%  Methanol /80%  Water 
PECLET  RF  PULSE 

5.0  1.0  0.322 

110 
0.440  0.0000 

0.491  0.0000 

0.542  0.0000 

0.593  0.0000 

0.645  0.0000 

0.696  0.0000 

0.747  0.0000 

0.798  0.0000 

0.850  0.0000 

0.901  0.0000 

0.952  0.0000 

1.003  0.0000 

1.055  0.0017 

1.106  0.0017 

1.157  0.0043 

1.208  0.0173 

1.260  0.0605 

1.311  0.1710 

1.362  0.3301 

1.413  0.5278 

1.465  0.7038 

1.516  0.7881 

1.567  0.8223 

1.618  0.7540 

1.670  0.7254 

1.721  0.5936 

1.772  0.4118 

1.823  0.1756 

1.875  0.0421 

1.926  0.0000 

1.977  0.0000 

2.028  0.0000 
2.080  0.0000 
2 . 131  •  0 . 0000 
2.182  0.0000 
2.233  0.0000 

2.285  0.0000 

2.336  0.0000 

2.387  0.0000 

2.438  0.0000 

C$stop 


68 


1 

2.1  25  40 

Methanol  Run:  50%  Methanol/50%  Water 
PECLET  RF  PULSE 


o 

in 

1. 

1  1 

0 

0.406 

0.0000 

0.458 

0.0000 

0.509 

0.0000 

0.560 

0.0000 

0.611 

0.0000 

0.663 

0.0000 

0.714 

0.0000 

0.765 

0.0000 

0.816 

0.0000 

0.868 

0.0000 

0.919 

0.0000 

0.970 

0.0000 

1.021 

0.0000 

1.073 

0.0000 

1.124 

0.0000 

1.175 

0.0177 

1.226 

0.0789 

1.278 

0.2270 

1.329 

0.3950 

1.380 

0.6083 

1.431 

0.7144 

1.483 

0.7184 

1.534 

0.7253 

1.585 

0.7146 

1.636 

0.6797 

1.688 

0.5980 

1.739 

0.4776 

1.790 

0.2658 

1.841 

0.0970 

1.893 

0.0184 

1.944 

0.0000 

1.995 

0.0000 

2.046 

0.0000 

2.098  - 

0.0000 

2.149 

0.0000 

2.200 

0.0000 

2.251 

0.0000 

2.303 

0.0000 

2.354 

0.0000 

2.405 

0.0000 

C$stop 


69 


1 

2  .  1 

25  40 

PCE  Run: 

0%  Methanol/100%  Water 

PECLET 

RF  PULSE 

159. 

1.0  0.334 

0  1 

0  0  0 

0.552 

0.0000 

0.603 

0.0000 

0.655 

0.0000 

0.706 

0.0000 

0.757 

0.0000 

0.808 

0.0000 

0.860 

0.0000 

0.911 

0.0000 

0.962 

0.0005 

1.013 

0.0021 

1.065 

0.0084 

1.116 

0.0597 

1.167 

0.2889 

1.218 

0.5989 

1.270 

0.8106 

1.321 

0.9129 

1.372 

0.9490 

1.423 

0.9235 

1.475 

0.7407 

1.526 

0.4700 

1.577 

0.2397 

1.628 

0.1175 

1.680 

0.0521 

1.731 

0.0261 

1.782 

0.0160 

1.833 

0.0123 

1.885 

0.0103 

1.936 

0.0089 

1.987 

0.0078 

2.038 

0.0067 

2.090 

0.0058 

2.141 

0.0048 

2.192 

0.0036 

2.243  - 

0.0028 

2.295 

0.0022 

2.346 

0.0017 

2.397 

0.0013 

2.448 

0.0011 

2.500 

0.0010 

2.551 

0.0008 

C$stop 


70 


1 


2  1 

25  40 

PCE  Run: 

1%  Methanol/99%  Water 

PECLET 

RF  PULSE 

153. 

1.0  0.337 

0  1 

0 

0.504 

0.0000 

0.555 

0.0001 

0.607 

0.0005 

0.658 

0.0013 

0.709 

0.0028 

0.760 

0.0057 

0.812 

0.0099 

0.863 

0.0165 

0 . 914 

0.0247 

0.965 

0.0359 

1.017 

0.0530 

1.068 

0.0820 

1 . 119 

0.1413 

1.170 

0.2641 

1.222 

0.4469 

1.273 

0.6510 

1.324 

0.8571 

1.375 

0.8900 

1.427 

0.8853 

1.478 

0.7760 

1.529 

0.6292 

1.580 

0.4164 

1.632 

0.2006 

1.683 

0.0725 

1.734 

0.0227 

1.785 

0.0105 

1.837 

0.0073 

1.888 

0.0059 

1.939 

0.0052 

1.990 

0.0050 

2.042 

0.0048 

2 . 093 

0.0047 

2.144 

0.0048 

2.195 

0.0050 

2.247 

0.0049 

2.298 

0.0048 

2.349 

0.0042 

2.400 

0.0038 

2.452 

0.0034 

2.503 

0.0031 

C$stop 


71 


1 


2  1 

25  40 

PCE  Run: 

5%  Methanol/95%  Water 

PECLET 

RF  PULSE 

305. 

1.0  0.309 

0  1 

0 

0 . 606 

0.0003 

0 . 658 

0.0007 

0.709 

0.0019 

0.760 

0.0041 

0.811 

0.0077 

0.863 

0.0133 

0.914 

0.0218 

0.965 

0.0338 

1 . 016 

0.0499 

1 . 068 

0.0785 

1 . 119 

0.1387 

1.170 

0.2247 

1.221 

0.3449 

1.273 

0.5036 

1.324 

0.7017 

1.375 

0.8171 

1.426 

0.8316 

1.478 

0.7858 

1.529 

0.6491 

1.580 

0.3957 

1.631 

0.1656 

1.683 

0.0522 

1.734 

0.0200 

1.785 

0.0109 

1.836 

0.0070 

1.888 

0.0055 

1.939 

0.0047 

1.990 

0.0041 

2.041 

0.0038 

2 . 093 

0.0036 

2 . 144 

0.0034 

2.195 

0.0032 

2 . 246 

0.0031 

2.298 

0.0029 

2.349 

'0.0028 

2.400 

0.0027 

2.451 

0.0026 

2.503 

0.0027 

2.554 

0.0025 

2.605 

0.0023 

C$stop 
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1 


2  1 

25  40 

PCE  Run: 

10%  Methanol/90%  Water 

PECLET 

RF  PULSE 

276. 

1.0  0.320 

0  1 

0 

0.450 

0.0002 

0.501 

0.0007 

0.553 

0.0018 

0.604 

0.0044 

0.655 

0.0077 

0.706 

0.0122 

0.758 

0.0181 

0.809 

0.0242 

0.860 

0.0328 

0.911 

0.0406 

0.963 

0.0507 

1.014 

0.0659 

1.065 

0.0934 

1.116 

0.1523 

1.168 

0.2284 

1.219 

0.3264 

1.270 

0.5001 

1.321 

0.6810 

1.373 

0.7971 

1.424 

0.8124 

1.475 

0.7504 

1.526 

0.6556 

1.578 

0.4617 

1.629 

0.2220 

1.680 

0.0719 

1.731 

0.0266 

1.783 

0.0139 

1.834 

0.0093 

1.885 

0.0068 

1.936 

0.0054 

1.988 

0.0045 

2.039 

0.0038 

2.090 

0.0033 

2.141 

0.0029 

2.193 

■  0.0026 

2.244 

0.0024 

2.295 

0.0022 

2.346 

0.0021 

2.398 

0.0020 

2.449 

0.0019 

C$stop 
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1 


2  1 

25  40 

PCE  Run; 

20%  Methanol/80%  Water 

PECLET 

RF  PULSE 

328. 

1.0  0.313 

0  1 

0 

0.483 

0.0125 

0.535 

0.0172 

0.586 

0.0221 

0.637 

0.0276 

0.688 

0.0330 

0.740 

0.0374 

0.791 

0.0444 

0.842 

0.0516 

0.893 

0.0607 

0.945 

0.0735 

0.996 

0.0930 

1.047 

0.1491 

1.098 

0.2285 

1.150 

0.3114 

1.201 

0.3575 

1.252 

0.4748 

1.303 

0.6591 

1.355 

0.8305 

1.406 

0.8706 

1.457 

0.7595 

1.508 

0.6463 

1.560 

0.4454 

1.611 

0.2223 

1.662 

0.0708 

1.713 

0.0333 

1.765 

0.0120 

1.816 

0.0140 

1.867 

0.0111 

1.918 

0.0092 

1.970 

0.0075 

2.021 

0.0072 

2.072 

0.0056 

2.123 

0.0052 

2.175 

0.0047 

2.226 

*0.0046 

2.277 

0.0044 

2.328 

0.0043 

2.380 

0.0044 

2.431 

0.0043 

2.482 

0.0040 

C$stop 


74 


1 


2  1 

25  40 

PCE  Run: 

50%  Methanol/50%  Water 

PECLET 

RF  PULSE 

221. 

1.0  0.326 

0  1 

0 

0.433 

0.0084 

0.485 

0.0136 

0.536 

0.0191 

0.587 

0.0237 

0.638 

0.0296 

0.690 

0.0352 

0.741 

0.0384 

0.792 

0.0453 

0.843 

0.0511 

0.895 

0.0757 

0.946 

0.0688 

0.997 

0.0902 

1.048 

0.1298 

1.100 

0.2016 

1.151 

0.2740 

1.202 

0.3354 

1.253 

0.4179 

1.305 

0.5865 

1.356 

0.8007 

1.407 

0.9018 

1.458 

0.8672 

1.510 

0.7320 

1.561 

0.4621 

1.612 

0.2397 

1.663 

0.1438 

1.715 

0.0839 

1.766 

0.0619 

1.817 

0.0458 

1.868 

0.0300 

1.920 

0.0191 

1.971 

0.0122 

2.022 

0.0083 

2.073 

0.0057 

2.125 

0.0048 

2.176 

■  0.0039 

2.227 

0.0034 

2.278 

0.0030 

2.330 

0.0026 

2.381 

0.0024 

2.432 

0.0022 

C$stop 
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Appendix  C.  CFITM  Output  Files 


76 


NON-LINEAR  LEAST  SQUARES  ANALYSIS 


EQUILIBRIUM  TRANSPORT  (MODEL  A) 
THIRD- TYPE  BOUNDARY  CONDITION 
Methanol  Run;  0%  Methanol/100%  Water 


INITIAL 

VALUES 

OF  COEFFICIENTS 

NO 

NAME 

INITIAL  VALUE 

1 

PECLET 

5.000 

2 

RF 

1.000 

3 

OBSERVED 

PUL 

1  DATA 

0.319 

OBS.  NO. 

PORE  VOLUME 

CONCENTRATION 

1 

0.4000 

0.0000 

2 

0.4510 

0.0000 

3 

0.5030 

0.0000 

4 

0.5540 

0.0000 

5 

0.6050 

0.0000 

S 

0.6560 

0.0000 

7 

0.7080 

0.0000 

8 

0.7590 

0.0000 

9 

0.8100 

0.0000 

10 

0.8610 

0.0000 

11 

0.9130 

0.0000 

12 

0.9640 

0.0028 

13 

1.0150 

0.0111 

14 

1.0660 

0.0184 

15 

1.1180 

0.0320 

16 

1.1690 

0.0563 

17 

1.2200 

0.0911 

18 

1.2710 

0.1527 

19 

1.3230 

0.2566 

20 

1.3740 

0.4205 

21 

1.4250 

0.5675 

22 

1.4760 

0.6260 

23 

1.5280 

0.6289 

24 

1.5790 

0,5736 

25 

1.6300 

0.4732 

26 

1.6810 

0.3275 

27 

1.7330 

0.1455 

28 

1.7840 

0.0284 

29 

1.8350 

0.0000 

30 

1.8860 

0.0000 

31 

1.9380 

0.0000 

32 

1.9890 

0.0000 

33 

2.0400 

0.0000 

34 

2.0910 

0.0000 

35 

2.1430 

0.0000 

36 

“ 

2.1940 

0.0000 

37 

2.2450 

0.0000 

38 

2.2960 

0.0000 

77 


39 

40 


2.3480  0.0000 

2.3990  0.0000 


1 

2 


ITERATION 

SSQ 

PECLET 

0 

1.7375172 

5.00000 

1 

1.3608561 

3.95501 

2 

1.1148545 

10.48430 

3 

0.9179355 

21.69042 

4 

0.6162016 

39.98299 

5 

0.4030838 

66.54190 

6 

0.2705094 

100.95358 

7 

0.2326818 

130.64368 

8 

0.2254175 

148.79398 

9 

0.2245778 

156.06346 

10 

0.2245075 

158.40369 

11 

0.2245041 

159.08261 

12 

0.2245041 

159.09472 

CORRELATION 

MATRIX 

1 

2 

1  1.0000 

2  0.0432  1.0000 

NON-LINEAR 

LEAST  SQUARES 

ANALYSIS,  F: 

\BLE  NAME 

VALUE 

S.E.COEFF. 

PECLET 

159.09472 

23.1279 

RF 

1.34607 

0.0110 

RF 

1.00000 

1.61202 

1.57499 

1.26549 

1.41942 

1.31392 

1.35637 

1.34301 

1.34595 

1.34598 

1.34605 

1.34607 

1.34607 


95%  CONFIDENCE  LIMITS 
T-VALUE  LOWER  UPPER 

6.88  112.2719  205.9175 

122.14  1.3238  1.3684 


ORDERED  BY  COMPUTER  INPUT 


PORE  CONCENTRATION  RESI- 

NO  VOLUME  OBS.  FITTED 

DUAL 

PORE 

NO 

1 

0.400 

0.000 

0.000 

0.000 

13 

2 

0.451 

0.000 

0.000 

0.000 

12 

3 

0,503 

0.000 

0.000 

0.000 

14 

4 

0.554 

0.000 

0.000 

0.000 

1 

5 

0.605 

0.000 

0.000 

0.000 

2 

6 

0.656 

0.000 

0.000 

0.000 

3 

7 

0.708 

0.000 

0.000 

0.000 

4 

8 

0.759 

0.000 

0.000 

0.000 

5 

9 

0.810 

0.000 

0.000 

0.000 

6 

10 

0.861 

0.000 

0.000 

0.000 

7 

11 

0.913 

0.000 

0.000 

0.000 

8 

12 

0.964 

0.003 

0.001 

0.001 

9 

13 

1.015 

0.011 

0.006 

0.005 

10 

14 

1.066 

0.018 

0.018 

0.000 

40 

15 

1.118 

0.032 

0.048 

-0.016 

39 

16 

1.169 

0.056 

0.103 

-0.047 

11 

17 

1.220 

0.091 

0.189 

-0.098 

38 

18 

1.271 

0.153 

0.303 

-0.150 

37 

19 

1.323 

0.257 

0.434 

-0.178 

36 

20 

1.374 

0.420 

‘  0.558 

-0.138 

35 

21 

1.425 

0.568 

0.656 

-0.088 

34 

22 

1.476 

0.626 

0.708 

-0.082 

33 

. ORDERED  BY  RESIDUALS- - 

CONCENTRATION  RES I - 

VOLUME  OBS.  FITTED  DUAL 


1.015 

0.011 

0.964 

0.003 

1.066 

0.018 

0.400 

0.000 

0.451 

0,000 

0.503 

0.000 

0.554 

0.000 

0.605 

0.000 

0.656 

0.000 

0.708 

0.000 

0.759 

0.000 

0.810 

0.000 

0.861 

0.000 

2.399 

0.000 

2.348 

0.000 

0,913 

0.000 

2.296 

0.000 

2.245 

0.000 

2.194 

0.000 

2.143 

0.000 

2.091 

0.000 

2.040 

0.000 

0.006 

0.005 

0.001 

0.001 

0.018 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.001 

-0.001 

0.001 

-0.001 

0.003 

-0.003 

0.007 

-0.007 

0.014 

-0.014 

78 


23 

1.528 

0.629 

0.704 

-0.075 

24 

1.579 

0.574 

0.646 

-0.073 

25 

1.630 

0.473 

0.550 

-0.077 

26 

1.681  • 

0,328 

0.435 

-0.108 

27 

1.733 

0.145 

0.318 

-0.173 

28 

1.784 

0.028 

0.219 

-0.190 

29 

1.835 

0.000 

0.141 

-0.141 

30 

1.886 

0.000 

0.086 

-0.086 

31 

1.938 

0.000 

0.049 

-0.049 

32 

1.989 

0.000 

0.026 

-0.026 

33 

2.040 

0.000 

0.014 

-0.014 

34 

2.091 

0.000 

0.007 

-0.007 

35 

2.143 

0.000 

0.003 

-0.003 

36 

2.194 

0.000 

0.001 

-0.001 

37 

2.245 

0.000 

0.001 

-0.001 

38 

2.296 

0.000 

0.000 

0.000 

39 

2.348 

0.000 

0.000 

0.000 

40 

2.399 

0.000 

0.000 

0.000 

15 

1.118 

0.032 

0.048 

-0.016 

32 

1.989 

0.000 

0.026 

-0.026 

16 

1.169 

0.056 

0.103 

-0.047 

31 

1.938 

0.000 

0.049 

-0.049 

24 

1.579 

0.574 

0.646 

-0.073 

23 

1.528 

0.629 

0.704 

-0.075 

25 

1.630 

0.473 

0.550 

-0.077 

22 

1.476 

0.626 

0.708 

-0.082 

30 

1.886 

0.000 

0.066 

-0.086 

21 

1.425 

0.568 

0.656 

-0.088 

17 

1.220 

0.091 

0.189 

-0.098 

26 

1.681 

0.328 

0.435 

-0.108 

20 

1.374 

0.420 

0.558 

-0.138 

29 

1.835 

0.000 

0.141 

-0.141 

18 

1,271 

0.153 

0.303 

-0.150 

27 

1.733 

0.145 

0,318 

-0.173 

19 

1.323 

0.257 

0.434 

-0.178 

28 

1.784 

0.028 

0.219 

-0.190 

END  OF  PROBLEM 


NON-LINEAR  LEAST  SQUARES  ANALYSIS 

EQUILIBRIUM  TRANSPORT  (MODEL  A) 
THIRD-TYPE  BOUNDARY  CONDITION 
Methanol  Run;  IV  Methanol/99V  Water 


INITIAL  VALUES  OF  COEFFICIENTS 


NO 

1 

2 

3 


NAME 

PECLET 

RF 

PUL 


INITIAL  VALUE 
5.000 
1.000 
0.326 


OBSERVED  DATA 


OBS.  NO. 
1 
2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 
17 
16 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 


PORE  VOLUME 
0.4310 
0.4830 
0.5340 
0.5850 
0.6360 
0.6880 
0.7390 
0,7900 
0.8410 
0.8930 
0.9440 
0.9950 
1.0460 
1.0980 
1.1490 
1.2000 
1.2510 
1.3030 
1.3540 
1.4050 
1.4560 
1.5080 
1.5590 
1.6100 
1.6610 
1.7130 
1.7640 
1.8150 
1.8660 
1.9180 
1.9690 
2.0200 
2.0710 
2.1230 
2.1740 
2.2250 
2.2760 
2.3280 


CONCENTRATION 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0046 

0.0136 

0.0245 

0.0436 

0.0620 

0.1352 

0.2179 

0.3383 

0.4895 

0.6194 

0.6339 

0.6113 

0.5364 

0.4029 

0.2162 

0.0678 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 
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1 


39 

2.3790 

40 

2.4300 

ITERATION 

SSQ 

0 

1.7593414 

1 

1.3673960 

2 

1.1304704 

3 

0.9330030 

4 

0.6351175 

5 

0.4244203 

6 

0.2966283 

7 

0.2591371 

8 

0.2515045 

9 

0.2505274 

10 

0.2S043S1 

11 

0.2504297 

12 

0.2504297 

CORRELATION 

MATRIX 

1 

2 

1  1.0000 

2  0.0425  1.0000 

0.0000 

0.0000 


PECLET 

RF 

5.00000 

1.00000 

4.02193 

1.60697 

10.36198 

1.57879 

20.88302 

1.26978 

38.60093 

1.41651 

63.76383 

1.31535 

96.48217 

1.35419 

124.63864 

1.34265 

142.14987 

1.34461 

149.52063 

1.34467 

152.03288 

1.34476 

152.80571 

1.34478 

152.83439 

1.34478 

NON-LINEAR  LEAST  SQUARES  ANALYSIS,  FINAL  RESULTS 


VARIABLE  NAME 

1  PECLET 

2  RF 


VALUE 

152.83439 

1.34478 


S . E . COEFF . 
23.2397 
0.0117 


95V  CONFIDENCE  LIMITS 
T- VALUE  LOWER  UPPER 

6.58  105.7853  199.8834 

114.60  1.3210  1.3685 


. ORDERED  BY  COMPUTER 

PORE  CONCENTRATION  RES I- 

INPUT-- 

PORE 

NO 

VOLUME 

OBS. 

FITTED 

DUAL 

NO 

1 

0.431 

0.000 

0.000 

0.000 

12 

2 

0.483 

0.000 

0.000 

0.000 

13 

3 

0.534 

0.000 

0.000 

0.000 

1 

4 

0.585 

0.000 

0.000 

0.000 

2 

5 

0.636 

0.000 

0.000 

0.000 

3 

6 

0.688 

0.000 

0.000 

0.000 

4 

7 

0.739 

0.000 

0.000 

0.000 

5 

8 

0.790 

0.000 

0.000 

0.000 

6 

9 

0.841 

0.000 

0.000 

0.000 

7 

10 

0.893 

0.000 

0.000 

0.000 

8 

11 

0.944 

0.000 

0.001 

-0.001 

9 

12 

0.995 

0.005 

0.004 

0.001 

40 

13 

1.046 

0.014 

0.014 

0.000 

39 

14 

1.098 

0.025 

0.037 

-0.013 

10 

15 

1.149 

0.044 

0.084 

-0.040 

38 

16 

1.200 

0.082 

0.159 

-0.077 

37 

17 

1.251 

0.135 

0.262 

-0.127 

11 

18 

1.303 

0.218 

0.388 

-0.170 

36 

19 

1.354 

0.3.38 

0.515 

-0.176 

35 

20 

1.405 

0.489 

0.623 

-0.133 

34 

21 

1.456 

0.619 

0.694 

-0.074 

33 

22 

1.508 

0.634 

0.714 

-0.080 

14 

. ORDERED  BY  RESIDUALS- - 

CONCENTRATION  RES I - 
VOLUME  OBS .  FITTED  DUAL 


0.995 

0.005 

1.046 

0.014 

0.431 

0.000 

0.483 

0.000 

0.534 

0.000 

0.585 

0.000 

0.636 

0.000 

0.688 

0.000 

0.739 

0.000 

0.790 

0.000 

0.841 

0.000 

2.430 

0.000 

2.379 

0.000 

0.893 

0.000 

2.328 

0.000 

2.276 

0.000 

0.944 

0.000 

2.225 

0.000 

2.174 

0.000 

2.123 

0.000 

2.071 

0.000 

1.098 

0.025 

0.004 

0.001 

0.014 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.001 

-0.001 

0.001 

-0.001 

0.001 

-0.001 

0.003 

-0.003 

0.005 

-0.005 

0.011 

-0.011 

0.037 

-0.013 
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23 

1.559 

0.611 

24 

1.610 

0.536 

25 

1.661 

0.403 

26 

1.713  . 

0.216 

27 

1.764 

0.068 

28 

1.815 

0.000 

29 

1.866 

0.000 

30 

1.918 

0.000 

31 

1.969 

0.000 

32 

2.020 

0.000 

33 

2.071 

0.000 

34 

2.123 

0.000 

35 

2.174 

0.000 

36 

2.225 

0.000 

37 

2.276 

0,000 

38 

2.328 

0.000 

39 

2.379 

0.000 

40 

2.430 

0.000 

0.660 

-0.068 

32 

0.601 

-0.064 

31 

0.494 

-0.091 

15 

0.377 

-0.160 

24 

0.270 

-0.202 

23 

0.182 

-0.182 

30 

0.115 

-0.115 

21 

0.068 

-0.066 

16 

0.039 

-0.039 

22 

0.021 

-0.021 

25 

0.011 

-0.011 

29 

0.005 

-0.005 

17 

0.003 

-0.003 

20 

0.001 

-0.001 

26 

0.001 

-0.001 

18 

0.000 

0.000 

19 

0.000 

0.000 

28 

0.000 

0.000 

27 

2.020 

0.000 

0.021 

1.969 

0.000 

0.039 

1.149 

0.044 

0.084 

1.610 

0.536 

0.601 

1.559 

0.611 

0.680 

1.918 

0.000 

0.068 

1.456 

0.619 

0.694 

1.200 

0.082 

0.159 

1.508 

0.634 

0.714 

1.661 

0.403 

0.494 

1.866 

0.000 

0.115 

1.251 

0.135 

0.262 

1.405 

0.489 

0.623 

1.713 

0.216 

0.377 

1.303 

0.218 

0.388 

1.354 

0.338 

0.515 

1.815 

0.000 

0.182 

1.764 

0.068 

0.270 

END  OF  PROBLEM 


0.021 

0.039 

0.040 

0.064 

0.066 

0.068 

0.074 

0.077 

0.080 

0.091 

0.115 

0.127 

0.133 

0.160 

0.170 

0.176 

0.182 

0.202 
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NON-LINEAR  LEAST  SQUARES  ANALYSIS 


EQUILIBRIUM  TRANSPORT  (MODEL  A) 
THIRD-TYPE  BOUNDARY  CONDITION 
Methanol  Run;  5V  Methanol/95%  Water 


INITIAL  VALUES  OF  COEFFICIENTS 


0  NAME  INITIAL  VALUE 

1  PECLET  5.000 

2  RF  1,000 

3  PUL  0.326 


OBSERVED  DATA 


!.  NO. 

PORE  VOLUME 

1 

0.4560 

2 

0.5080 

3 

0.5590 

4 

0.6100 

5 

0.6610 

6 

0.7130 

7 

0.7640 

8 

0.8150 

9 

0.6660 

10 

0.9180 

11 

0.9690 

12 

1.0200 

13 

1.0710 

14 

1.1230 

15 

1.1740 

16 

1.2250 

17 

1.2760 

18 

1.3280 

19 

1.3790 

20 

1.4300 

21 

1.4810 

22 

1.5330 

23 

1.5840 

24 

1.6350 

25 

1.6860 

26 

1.7380 

27 

1.7890 

28 

1.8400 

29 

1.8910 

30 

1.9430 

31 

1.9940 

32 

2.0450 

33 

2.0960 

34 

2.1480 

35 

2.1990 

36 

2.2500 

37 

2.3010 

38 

2.3530 

CONCENTRATION 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0024 

0.0108 

0.0318 

0.0926 

0.2447 

0.4720 

0.6573 

0.7483 

0.7455 

0.7510 

0.6581 

0.4568 

0.1993 

0.0302 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 
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39 

40 


2.4040 

2.4550 


0.0000 

0.0000 


1 


ITERATION  SSQ 

0  2.5125510 

1  2.0375895 

2  1.8661717 

3  1.3585833 

4  1.2421442 

5  0.6307918 

6  0.3239347 

7  0.1571142 

8  0.1274750 

9  0.1242359 

10  0.1240597 

11  0.1240537 

12  0.1240537 

13  0.1240537 
CORRELATION  MATRIX 


1  2 

1  1.0000 

2  0.0775  1.0000 


PECLET 

RF 

5.00000 

1.00000 

6.02197 

1.63007 

16.81566 

1.14096 

28.57618 

1.46650 

55.70067 

1.18348 

77.18728 

1.37174 

132.13273 

1.28056 

206.50526 

1.32136 

264.84063 

1.31176 

294.52444 

1.31329 

302.77202 

1.31336 

304.56661 

1.31339 

304.74452 

1.31339 

304.75780 

1.31339 

NON-LINEAR  LEAST  SQUARES  ANALYSIS,  FINAL  RESULTS 


VARIABLE  NAME  VALUE  S . E . COEFF . 

1  PECLET  304.75780  39.7658 

2  RF  1.31339  0.0059 


95V  CONFIDENCE  LIMITS 
T-VALUE  LOWER  UPPER 

7.66  224.2514  305.2642 

220.81  1.3014  1.3254 


. ORDERED  BY  COMPUTER  INPUT 

PORE  CONCENTRATION  RESI- 
NO  VOLUME  OBS .  FITTED  DUAL 

1  0.456  0.000  O.'OOO  0.000 

2  0.508  0.000  0.000  0.000 

3  0.559  0.000  0.000  0.000 

4  0.610  0.000  0.000  0.000 

5  0.661  0.000  0.000  0.000 

6  0.713  0.000  0.000  0.000 

7  0.764  0.000  0.000  0.000 

8  0.815  0.000  0.000  0.000 

9  0.866  0.000  0.000  0.000 

10  0.918  0.000  0.000  0.000 

11  0.969  0.000  0.000  0.000 

12  1.020  0.000  0.001  -0.001 

13  1.071  0.002  0.006  -0.003 

14  1.123  0.011  0.026  -0.015 

15  1.174  0,032  0.083  -0.051 

16  1.225  0.093  0.194  -0.102 

17  1.276  0.245  0.360  -0.116 

18  1.328  0-.472  0.554  -0.082 

19  1.379  0.657  0.724  -0.066 

20  1.430  0.748  0.838  -0.090 

21  1.481  0.746  0.875  -0.130 


- ORDERED  BY  RESIDUALS - 

PORE  CONCENTRATION  RES  I  - 

NO  VOLUME  OBS.  FITTED  DUAL 
1  0.456  0,000  0.000  0.000 

2  0.508  0.000  0.000  0.000 

3  0.559  0.000  0.000  0.000 

4  0.610  0.000  0.000  0.000 

5  0.661  0.000  0.000  0.000 

6  0.713  0.000  0.000  0.000 

7  0.764  0.000  0.000  0.000 

40  2.455  0.000  0.000  0.000 

8  0.815  0.000  0.000  0.000 

39  2.404  0.000  0.000  0.000 

38  2.353  0.000  0.000  0.000 

9  0.866  0.000  0.000  0.000 

37  2.301  0.000  0.000  0.000 

36  2.250  0.000  0.000  0.000 

10  0.918  0.000  0.000  0.000 

35  2.199  0.000  0.000  0.000 

34  2.148  0.000  0.000  0.000 

11  0.969  0.000  0.000  0.000 

33  2.096  0.000  0.000  0.000 

32  2.045  0.000  0.000  0.000 

12  1.020  0.000  0.001  -0.001 
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22 

1.533 

0.751 

23 

1.584 

0.658 

24 

1.635 

0.457 

25 

1.686 

0.199 

26 

1.738 

0.030 

27 

1.789 

0.000 

28 

1.840 

0.000 

29 

1.891 

0.000 

30 

1.943 

0.000 

31 

1.994 

0.000 

32 

2.045 

0.000 

33 

2.096 

0.000 

34 

2.148 

0.000 

35 

2.199 

0.000 

36 

2.250 

0.000 

37 

2.301 

0.000 

38 

2.353 

0.000 

39 

2.404 

0.000 

40 

2.455 

0.000 

0.624 

-0.073 

31 

0.693 

-0.035 

13 

0.513 

-0.056 

30 

0.332 

-0.133 

29 

0.185 

-0.155 

14 

0,091 

-0.091 

23 

0,039 

-0.039 

28 

0.015 

-0,015 

15 

0.005 

-0.005 

24 

0.002 

-0.002 

19 

0.000 

0.000 

22 

0.000 

0.000 

18 

0.000 

0.000 

20 

0.000 

0.000 

27 

0.000 

0.000 

16 

0.000 

0.000 

17 

0.000 

0.000 

21 

0,000 

0.000 

25 

0.000 

0.000 

26 

1.994 

0.000 

0.002 

1.071 

0.002 

0.006 

1.943 

0.000 

0.005 

1.891 

0.000 

0.015 

1.123 

0.011 

0.026 

1.584 

0.658 

0.693 

1.840 

0,000 

0.039 

1.174 

0.032 

0.083 

1.635 

0.457 

0.513 

1.379 

0.657 

0.724 

1.533 

0.751 

0.824 

1.328 

0.472 

0.554 

1.430 

0.748 

0.638 

1,789 

0.000 

0.091 

1.225 

0.093 

0.194 

1.276 

0.245 

0.360 

1.481 

0.746 

0.875 

1.686 

0.199 

0.332 

1.738 

0.030 

0.185 

END  OF  PROBLEM 


0.002 

0.003 

0.005 

0.015 

0.015 

0.035 

0.039 

0.051 

0.056 

0.066 

0.073 

0.082 

0.090 

0.091 

0.102 

0.116 

0.130 

0.133 

0.155 
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NON-LINEAR  LEAST  SQUARES  ANALYSIS 


EQUILIBRIUM  TRANSPORT  (MODEL  A) 
THIRD-TYPE  BOUNDARY  CONDITION 
Methanol  Run:  10%  Methanol/90%  Water 


INITIAL  VALUES  OF  COEFFICIENTS 


NO 

NAME 

INITIAL  VALUE 

1 

PECLET 

5.000 

2 

RF 

1.000 

3 

PUL 

0,328 

OBSERVED 

DATA 

OBS.  NO. 

PORE  VOLUME 

CONCENTRATION 

1 

0.3940 

0.0000 

2 

0.4450 

0.0000 

3 

0.4960 

0.0000 

4 

0.5480 

0.0000 

5 

0.5990 

0,0000 

6 

0.6500 

0.0000 

7 

0.7010 

0.0000 

8 

0.7530 

0.0000 

9 

0.8040 

O.OOOQ 

10 

0.8550 

0.0000 

11 

0.9060 

0.0000 

12 

0.9580 

0.0000 

13 

1.0090 

0,0000 

14 

1.0600 

0.0018 

15 

1.1110 

0.0041 

16 

1.1630 

0.0353 

17 

1.2140 

0.1728 

18 

1.2650 

0.4012 

19 

1.3160 

0.6129 

20 

1.3680 

0.7223 

21 

1.4190 

0.8058 

22 

1.4700 

0.7435 

23 

1.5210 

0.6161 

24 

1.5730 

0.5071 

25 

1.6240 

0.3088 

26 

1.6750 

0.0880 

27 

1.7260 

0.0000 

28 

1.7780 

0.0000 

29 

1.8290 

0.0000 

30 

1.8800 

0.0000 

31 

1.9310 

0.0000 

32 

1.9830 

0.0000 

33 

2.0340 

0.0000 

34 

2.0850 

0.0000 

35 

2.1360 

0.0000 

36 

2.1880 

0.0000 

37 

2.2390 

0.0000 

38 

2.2900 

0.0000 

1 


39 

2.3410 

40 

2.3930 

ITERATION 

SSQ 

0 

2.3998257 

1 

1.9379348 

2 

1.88SB709 

3 

X.32S7426 

4 

1.1300019 

5 

0.5767929 

6 

0.2952515 

7 

0.1622297 

8 

0.1413681 

9 

0.1394053 

10 

0.1393203 

11 

0.1393184 

12 

0.1393184 

CORRELATION 

MATRIX 

1 

2 

1  1.0000 

2  0.0764 

t  1.0000 

0.0000 

0.0000 


PECLET 

RF 

5.00000 

1.00000 

6.79970 

1.57860 

17.65172 

1.07969 

25.74722 

1.40975 

52.26530 

1.15333 

80.24640 

1.32332 

129.93195 

1.23941 

199.68445 

1.27236 

247.17486 

1.26468 

269.29171 

1.26545 

274.57353 

1.26543 

275.58952 

1.26542 

275.61210 

1.26542 

NON-LINEAR  LEAST  SQUARES  ANALYSIS,  FINAL  RESULTS 


VARIABLE  NAME  VALUE 

1  PECLET  275.61210 

2  RF  1.26542 


S.E.COEFF. 

37.8166 

0.0064 


95%  CONFIDENCE  LIMITS 
T-VALUE  LOWER  UPPER 

7.29  199.0518  352.1724 

199.14  1.2526  1.2783 


. ORDERED  BY  COMPUTER  INPUT . 

PORE  CONCENTRATION  RESI-  PORE 

NO  VOLUME  OBS .  FITTED  DUAL  NO 


1 

0.394 

0.000 

2 

0.445 

0.000 

3 

0.496 

0.000 

4 

0.548 

0.000 

5 

0.599 

0.000 

6 

0.650 

0.000 

7 

0.701 

0.000 

8 

0.753 

0.000 

9 

0.804 

0.000 

10 

0.855 

0.000 

11 

0.906 

0.000 

12 

0.958 

0.000 

13 

1.009 

0.000 

14 

1.060 

0.002 

15 

1.111 

0.004 

16 

1.163 

0.035 

17 

1.214 

0.173 

18 

1.265 

0.401 

19 

1.316 

0.613 

20 

1.368 

0.722 

21 

1.419 

0.806 

22 

1.470 

0.744 

0.000 

0.000 

1 

0.000 

0.000 

2 

0.000 

0.000 

3 

0.000 

0.000 

4 

0.000 

0.000 

5 

0.000 

0.000 

6 

0.000 

0.000 

7 

0.000 

0.000 

8 

0.000 

0.000 

40 

0.000 

0.000 

39 

0.000 

0.000 

9 

0.001 

-0.001 

38 

0.004 

-0.004 

37 

0.018 

-0.017 

10 

0.063 

-0.059 

36 

0.160 

-0.125 

35 

0.313 

-0.140 

11 

0.498 

-0.097 

34 

0.676 

-0.063 

33 

0.810 

-0.088 

12 

0.871 

-0.065 

32 

0.848 

-0.104 

31 

. ORDERED  BY  RESIDUALS- - 

CONCENTRATION  RES I - 
VOLUME  OBS .  FITTED  DUAL 


0.394 

0.000 

0.445 

0.000 

0.496 

0.000 

0.548 

0.000 

0.599 

0.000 

0.650 

0.000 

0.701 

0.000 

0.753 

0.000 

2.393 

0.000 

2.341 

0.000 

0.804 

0.000 

2.290 

0.000 

2.239 

0.000 

0.855 

0.000 

2,188 

0.000 

2.136 

0.000 

0.906 

0.000 

2.085 

0.000 

2.034 

0.000 

0.958 

0.000 

1.983 

0.000 

1.931 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0,000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0,000 

0.000 

0.000 

0.000 

0,000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.001 

-0.001 

0.001 

-0.001 

0.003 

-0,003 
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23 

1.521 

0.616 

24 

1.573 

0.507 

25 

1.624 

0.309 

26 

1.675* 

0.088 

27 

1.726 

0.000 

28 

1.778 

0.000 

29 

1.829 

0.000 

30 

1.880 

0.000 

31 

1.931 

0.000 

32 

1.983 

0.000 

33 

2.034 

0.000 

34 

2.085 

0.000 

35 

2.136 

0.000 

36 

2.188 

0.000 

37 

2.239 

0.000 

38 

2.290 

0.000 

39 

2.341 

0.000 

40 

2.393 

0.000 

0.741  -0.125  13 
0.571  -0.064  30 
0.388  -0.079  14 
0.231  -0,143  29 
0.120  -0.120  28 
0.055  -0.055  15 
0.022  -0.022  19 
0.008  -0.008  24 
0.003  -0.003  21 
0.001  -0.001  25 
0.000  0.000  20 
0.000  0.000  18 
0.000  0.000  22 
0.000  0.000  27 
0.000  0.000  23 
0.000  0.000  16 
0.000  0.000  17 
0.000  0.000  26 


1.009 

0.000 

0.004 

1.880 

0.000 

0.008 

1.060 

0.002 

0.018 

1.829 

0.000 

0.022 

1.778 

0.000 

0.055 

1.111 

0.004 

0.063 

1.316 

0.613 

0.676 

1.573 

0.507 

0.571 

1.419 

0.806 

0.871 

1.624 

0.309 

0.388 

1.368 

0.722 

0.810 

1.265 

0.401 

0.498 

1.470 

0.744 

0.848 

1.726 

0.000 

0.120 

1.521 

0.616 

0.741 

1.163 

0.035 

0.160 

1.214 

0.173 

0.313 

1.675 

0.088 

0.231 

END  OF  PROBLEM 


0.004 

0.008 

0.017 

0.022 

0.055 

0.059 

0.063 

0.064 

0.065 

0.079 

0.088 

0.097 

0.104 

0.120 

0.125 

0.125 

0.140 

>0.143 
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NON-LINEAR  LEAST  SQUARES  ANALYSIS 

EQUILIBRIUM  TRANSPORT  (MODEL  A) 
THIRD-TYPE  BOUNDARY  CONDITION 
Methanol  Run:  20V  Methanol/8 OV  Water 


INITIAL  VALUES  OF  COEFFICIENTS 


NO  NAME 

1  PECLET 

2  RF 

3  PUL 


INITIAL  VALUE 
5.000 
1.000 
0.322 


OBSERVED  DATA 


;.  NO. 

PORE  VOLUME 

1 

0.4400 

2 

0.4910 

3 

0.5420 

4 

0.5930 

5 

0.6450 

6 

0.6960 

7 

0.7470 

8 

0.7980 

9 

0.8500 

10 

0.9010 

11 

0.9520 

12 

1.0030 

13 

1.0550 

14 

1,1060 

15 

1.1570 

16 

1.2080 

17 

1.2600 

18 

1.3110 

19 

1.3620 

20 

1.4130 

21 

1.4650 

22 

1.5160 

23 

1.5670 

24 

1.6180 

25 

1.6700 

26 

1.7210 

27 

1.7720 

28 

1.8230 

29 

1.8750 

30 

1.9260 

31 

1.9770 

32 

2.0280 

33 

2.0800 

34 

2.1310 

35 

2.1820 

36 

2.2330 

37 

2.2850 

38 

2.3360 

CONCENTRATION 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0017 

0.0017 

0.0043 

0.0173 

0.0605 

0.1710 

0.3301 

0.5278 

0.7038 

0.7881 

0.8223 

0.7540 

0.7254 

0.5936 

0.4118 

0.1756 

0.0421 

0.0000 

0.0000 

0.0000 

0,0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 
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39 

2.3870 

0.0000 

40 

2.4380 

0.0000 

ITERATION 

SSQ 

PECLET 

RF 

0 

3.1821863 

5.00000 

1 . 00000 

1 

2.6683041 

4.62117 

1.76768 

2 

2.1786787 

11.35423 

1.40497 

3 

1.2182719 

37.32405 

1.44531 

4 

0.4977078 

103.62416 

1.35778 

5 

0.1255254 

198.04523 

1.43351 

6 

0.0311304 

282.26962 

1.40335 

7 

0.0222684 

323.83838 

1.41066 

8 

0.0222296 

328.22786 

1.41062 

9 

0.0222295 

328.24006 

1.41065 

CORRELATION  MATRIX 


1  2 

1  1.0000 

2  0.0764  1.0000 


NON-LINEAR  LEAST  SQUARES  ANALYSIS,  FINAL  RESULTS 


VARIABLE  NAME 

1  PECLET 

2  RF 


VALUE 

328.24006 

1.41065 


S.E.COEFF. 

17.6087 

0.0026 


95V  CONFIDENCE  LIMITS 
T-VALUE  LOWER  UPPER 

10.64  292.5910  363.8891 

543.86  1.4054  1.4159 


ORDERED  BY  COMPUTER  INPUT' 


PORE  CONCENTRATION  RESI- 

NO  VOLUME  OBS .  FITTED 

DUAL 

PORE 

NO 

1 

0.440 

0.000 

0.000 

0.000 

26 

2 

0.491 

0.000 

0.000 

0.000 

27 

3 

0.542 

0.000 

0.000 

0.000 

21 

4 

0.593 

0.000 

0.000 

0.000 

25 

5 

0.645 

0.000 

0.000 

0.000 

20 

6 

0.696 

0.000 

0.000 

0.000 

19 

7 

0.747 

0.000 

0.000 

0.000 

13 

8 

0.798 

0.000 

0.000 

0.000 

14 

9 

0.850 

0.000 

0.000 

0.000 

1 

10 

0.901 

0.000 

0.000 

0.000 

2 

11 

0.952 

0.000 

0.000 

0.000 

3 

12 

1.003 

0.000 

0.000 

0.000 

4 

13 

1.055 

0.002 

0.000 

0.002 

5 

14 

1.106 

0.002 

0.001 

0.001 

6 

15 

1.157 

0.004 

0.005 

-0.001 

7 

16 

1.208 

0.017 

0.023 

-0.006 

8 

17 

1.260 

0.060 

0.074 

-0.013 

9 

18 

1.311 

0.171 

0.174 

-0.003 

10 

19 

1.362 

0.330 

0.326 

0.004 

40 

20 

1.413 

0.528 

0.508 

0.020 

11 

21 

1.465 

0.704 

0.683 

0.021 

39 

22 

1.516 

0-788 

0.806 

-0.018 

38 

23 

1.567 

0.822 

0.857 

-0.035 

12 

24 

1.618 

0.754 

0.823 

-0.069 

37 

25 

1.670 

0.725 

0.705 

0.020 

36 

. ORDERED  BY  RESIDUALS- - 

CONCENTRATION  RES I - 
VOLUME  OBS.  FITTED  DUAL 


1.721 

0.594 

1.772 

0.412 

1.465 

0.704 

1.670 

0.725 

1.413 

0.528 

1.362 

0.330 

1.055 

0.002 

1.106 

0.002 

0.440 

0.000 

0.491 

0.000 

0.542 

0.000 

0.593 

0.000 

0.645 

0.000 

0.696 

0.000 

0.747 

0.000 

0.798 

0.000 

0.850 

0.000 

0.901 

0.000 

2.438 

0.000 

0.952 

0.000 

2.387 

0.000 

2.336 

0.000 

1.003 

0.000 

2.285 

0.000 

2.233 

0.000 

0.537 

0.057 

0.360 

0.051 

0.683 

0,021 

0.705 

0.020 

0.508 

0.020 

0.326 

0.004 

0.000 

0.002 

0.001 

0.001 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 
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26 

1.721 

0.594 

27 

1.772 

0.412 

26 

1.823 

0.176 

29 

1.875 

0.042 

30 

1.926 

0.000 

31 

1.977 

0.000 

32 

2.028 

0.000 

33 

2.080 

0.000 

34 

2.131 

0.000 

35 

2.182 

0.000 

36 

2.233 

0.000 

37 

2.285 

0.000 

38 

2.336 

0.000 

39 

2.387 

0.000 

40 

2.438 

0.000 

0.537 

0.057 

35 

0.360 

0.051 

34 

0.212 

-0.037 

15 

0.109 

-0.066 

33 

0.050 

-0.050 

18 

0.020 

-0.020 

16 

0.007 

-0.007 

32 

0.002 

-0.002 

17 

0.001 

-0.001 

22 

0.000 

0.000 

31 

0.000 

0.000 

23 

0.000 

0.000 

28 

0.000 

0.000 

30 

0.000 

0.000 

29 

0.000 

0.000 

24 

2.182 

0.000 

0.000 

2.131 

0.000 

0.001 

1.157 

0.004 

0.005 

2.080 

0.000 

0.002 

1.311 

0.171 

0.174 

1.208 

0.017 

0.023 

2.028 

0.000 

0.007 

1.260 

0.060 

0.074 

1.516 

0.788 

0.806 

1.977 

0.000 

0.020 

1.567 

0.822 

0.857 

1.823 

0.176 

0.212 

1.926 

0.000 

0.050 

1.875 

0.042 

0.109 

1.618 

0.754 

0.823 

END  OF  PROBLEM 


0.000 

0.001 

0.001 

0.002 

0.003 

0.006 

0.007 

0.013 

0.018 

0.020 

0.035 

0.037 

0.050 

0.066 

0.069 
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NON-LINEAR  LEAST  SQUARES  ANALYSIS 


EQUILIBRIUM  TRANSPORT  (MODEL  A) 
THIRD-TYPE  BOUNDARY  CONDITION 
Methanol  Run:  50%  Methanol/50%  Water 


INITIAL  VALUES  OF  COEFFICIENTS 


fO  NAME  INITIAL  VALUE 

1  PECLET  5.000 

2  RF  1.000 

3  PUL  0.331 


OBSERVED  DATA 


OBS.  NO. 
1 
2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 
33 

36 

37 

38 


PORE  VOLUME 
0.4060 
0.4580 
0.5090 
0.5600 
0.6110 
0.6630 
0.7140 
0.7650 
0.8160 
0.8680 
0.9190 
0.9700 
1.0210 
1.0730 
1.1240 
1.1750 
1.2260 
1.2780 
1.3290 
1.3800 
1.4310 
1.4830 
1.5340 
1.5850 
1.6360 
1.6880 
1.7390 
1.7900 
1.8410 
1.8930 
1.9440 
1.9950 
2.0460 
2.0980 
2.1490 
2.2000 
2.2510 
2.3030 


CONCENTRATION 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0177 

0.0789 

0.2270 

0.3950 

0.6083 

0.7144 

0.7184 

0.7253 

0.7146 

0.6797 

0.5980 

0.4776 

0.2658 

0.0970 

0.0184 

0.0000 

0.0000 

0,0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 
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39 

2.3540 

0.0000 

40 

2.4050 

0.0000 

ITERATION 

SSQ 

PECLET 

RF 

0 

2.9583570 

5.00000 

1.00000 

1 

2.3833102 

5.80803 

1.73222 

2 

1.8883046 

13.34307 

1.30709 

3 

1,0679509 

39.74408 

1.47595 

4 

0.9037115 

91.46498 

1.25007 

5 

0.2358949 

114.46765 

1,40446 

6 

0.0712604 

182.02134 

1.35552 

7 

0.0467000 

218.05743 

1.37170 

8 

0.0465927 

221.29075 

1.37083 

9 

0.0465926 

221.15155 

1.37086 

10 

0.0465926 

221.15960 

1.37086 

1  CORRELATION  MATRIX 


1  2 

1  1.0000 

2  0.0615  1.0000 


NON-LINEAR  LEAST  SQUARES  ANALYSIS,  FINAL  RESULTS 


VARIABLE  NAME 

VALUE 

S.E.COEFF. 

T-VALUE 

95V  CONFIDENCE  LIMITS 
LOWER  UPPER 

1  PECLET 

221.15960 

15.5099 

14.26 

189.7597 

252.5595 

2  RF 

1.37086 

0.0043 

319.91 

1.3622 

1.3795 

ORDERED  BY  COMPUTER  INPUT . ORDERED  BY  RESIDUALS 


PORE  CONCENTRATION  RESI 

- 

PORE 

CONCENTRATION 

RESI 

- 

NO 

VOLUME 

OBS. 

FITTED 

DUAL 

NO 

VOLUME 

OBS. 

FITTED 

DUAL 

1 

0.406 

0.000 

0.000 

0.000 

27 

1.739 

0.478 

0.383 

0.094 

2 

0.458 

0.000 

0.000 

0.000 

20 

1.380 

0.606 

0.525 

0.083 

3 

0.509 

0.000 

0.000 

0.000 

26 

1.688 

0.598 

0.529 

0.069 

4 

0.560 

0.000 

0.000 

0.000 

21 

1.431 

0.714 

0.664 

0.050 

5 

0.611 

0.000 

0.000 

0.000 

19 

1.329 

0.395 

0,371 

0.024 

6 

0.663 

0.000 

0,000 

0.000 

25 

1.636 

0.680 

0.667 

0.013 

7 

0.714 

0.000 

0.000 

0.000 

26 

1.790 

0.266 

0.253 

0.012 

8 

0.765 

0.000 

0.000 

0.000 

1 

0.406 

0.000 

0.000 

0.000 

9 

0.816 

0.000 

0.000 

0.000 

2 

0.458 

0.000 

0.000 

0.000 

10 

0.868 

0.000 

0.000 

0.000 

3 

0.509 

0.000 

0.000 

0.000 

11 

0.919 

0.000 

0.000 

0.000 

4 

0.560 

0.000 

0.000 

0.000 

12 

0.970 

0.000 

0.000 

0.000 

5 

0.611 

0.000 

0.000 

0.000 

13 

1.021 

0.000 

0.001 

-0.001 

6 

0.663 

0.000 

0.000 

0.000 

14 

1.073 

0.000 

0.005 

-0.005 

7 

0.714 

0.000 

0.000 

0.000 

15 

1.124 

0.000 

0.018 

-0.018 

6 

0.765 

0.000 

0.000 

0.000 

16 

1.175 

0.018 

0.052 

-0.034 

9 

0.816 

0.000 

0.000 

0.000 

17 

1.226 

0.079 

0.119 

-0,041 

10 

0.666 

0.000 

0.000 

0.000 

18 

1.278 

0.227 

0.230 

-0.003 

40 

2.405 

0.000 

0.000 

0.000 

19 

1.329 

0.395 

0.371 

0.024 

11 

0.919 

0.000 

0.000 

0.000 

20 

1.380 

0.608 

0.525 

0.083 

39 

2.354 

0.000 

0.000 

0.000 

21 

1.431 

0^714 

0.664 

0.050 

38 

2.303 

0.000 

0.000 

0.000 

22 

1.483 

0.718 

0.763 

-0.045 

12 

0.970 

0.000 

0.000 

0.000 

23 

1.534 

0.725 

0.798 

-0.073 

37 

2.251 

0.000 

0.000 

0.000 

24 

1.585 

0.715 

0.763 

-0.049 

36 

2.200 

0.000 

0.001 

-0,001 
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25 

1.636 

0.680 

0.667 

0.013 

13 

1.021 

0.000 

0.001 

-0.001 

26 

1 .688 

0.598 

0.529 

0.069 

35 

2.149 

0.000 

0.001 

-0 . 001 

27 

1.739 

0.478 

0.383 

0.094 

18 

1.278 

0.227 

0.230 

-0 . 003 

28 

1.790 

0.266 

0.253 

0.012 

34 

2.098 

0.000 

0.004 

-0 . 004 

29 

1.841 

0.097 

0.153 

-0.056 

14 

1.073 

0.000 

0.005 

-0 . 005 

30 

1.893 

0.018 

0.084 

-0.066 

33 

2.046 

0.000 

0.009 

-0.009 

31 

1 . 944 

0.000 

0.043 

-0.043 

15 

1.124 

0.000 

0.018 

-0.018 

32 

1.995 

0.000 

0.020 

-0.020 

32 

1.995 

0.000 

0.020 

-0 . 020 

33 

2.046 

0.000 

0.009 

-0.009 

16 

1.175 

0.018 

0.052 

-0 . 034 

34 

2.098 

0.000 

0.004 

-0.004 

17 

1.226 

0.079 

0.119 

-0.041 

35 

2.149 

0.000 

0.001 

-0.001 

31 

1.944 

0.000 

0.043 

-0 . 043 

36 

2.200 

0.000 

0.001 

-0.001 

22 

1.483 

0.718 

0.763 

-0 . 045 

37 

2.251 

0.000 

0.000 

0.000 

24 

1.585 

0.715 

0.763 

-0.049 

38 

2.303 

0.000 

0.000 

0.000 

29 

1,841 

0.097 

0.153 

-0.056 

39 

2.354 

0.000 

0.000 

0.000 

30 

1.893 

0.018 

0.084 

-0 . 066 

40 

2,405 

0.000 

0.000 

0.000 

23 

1.534 

0.725 

0.798 

-0 . 073 

END  OF  PROBLEM 
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NON-LINEAR  LEAST  SQUARES  ANALYSIS 


EQUILIBRIUM  TRANSPORT  (MODEL  A) 
THIRD -TYPE  BOUNDARY  CX)NDITION 
PCE  Run:  0%  Methanol/100%  Water 


INITIAL  VALUES  OF  COEFFICIENTS 


NO  NAME 

1  PECLET 

2  RF 

3  PUL 


INITIAL  VALUE 
159.000 
1.000 
0.334 


OBSERVED  DATA 


OBS.  NO. 
1 
2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 
35. 

36 

37 

38 


PORE  VOLUME 
0.5520 
0.6030 
0.6550 
0.7060 
0.7570 
0.8080 
0.8600 
0.9110 
0.9620 
1.0130 
1.0650 
1.1160 
1.1670 
1.2180 
1.2700 
1.3210 
1.3720 
1.4230 
1.4750 
1.5260 
1.5770 
1.6280 
1.6800 
1.7310 
1.7820 
1.8330 
1.8850 
1.9360 
1.9870 
2.0380 
2.0900 
2.1410 
2.1920 
2.2430 
2.2950 
2.3460 
2.3970 
2.4480 


CONCENTRATION 

0.0000 

0,0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0005 

0.0021 

0.0084 

0.0597 

0.2889 

0.5989 

0.8106 

0.9129 

0.9490 

0.9235 

0.7407 

0.4700 

0.2397 

0.1175 

0.0521 

0.0261 

0.0160 

0.0123 

0.0103 

0.0089 

0.0078 

0.0067 

0.0058 

0.0048 

0.0036 

0.0028 

0.0022 

0.0017 

0.0013 

0.0011 
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39 

40 


2.5000 

2.5510 


0.0010 

0.0008 


ITERATION  SSQ  RF 

0  3.5356415  1.00000 

1  0.3227643  1.16641 

2  0.2397820  1.20674 

3  0.2355492  1.19223 

4  0.2327444  1.19765 

5  0.2327444  1.19765 

CORRELATION  MATRIX 


1 

1  1.0000 


NON-LINEAR  LEAST  SQUARES  ANALYSIS,  FINAL  RESULTS 


VARIABLE  NAME 
1  RF 


95V  CONFIDENCE  LIMITS 

VALUE  S.E.COEFF.  T- VALUE  LOWER  UPPER 

1.19765  0.0097  124.03  1.1781  1.2172 


ORDERED  BY  COMPUTER  INPUT-----  -----ORDERED  BY  RESIDUALS 


PORE  CONCENTRATION  RES I 

. 

PORE 

CONCENTRATION 

RES  I 

- 

NO 

VOLUME 

OBS. 

FITTED 

DUAL 

NO 

VOLUME 

OBS. 

FITTED 

DUAL 

1 

0.552 

0.000 

0.000 

0.000 

18 

1.423 

0.923 

0.741 

0.182 

2 

0.603 

0.000 

0.000 

0.000 

17 

1.372 

0.949 

0.788 

0.161 

3 

0.655 

0.000 

0.000 

0.000 

16 

1.321 

0.913 

0.768 

0.145 

4 

0.706 

0.000 

0.000 

0.000 

15 

1.270 

0.811 

0.686 

0.124 

5 

0.757 

0.000 

0.000 

0.000 

19 

1.475 

0.741 

0.637 

0.104 

6 

0.808 

0.000 

0.000 

0.000 

14 

1.218 

0.599 

0.557 

0.042 

7 

0.860 

0.000 

0.001 

-0.001 

30 

2.038 

0.007 

0.001 

0.006 

6 

0.911 

0.000 

0.007 

-0.007 

29 

1.987 

0.008 

0.002 

0.006 

9 

0.962 

0.001 

0.025 

-0.024 

31 

2.090 

0.006 

0.000 

0.006 

10 

1.013 

0.002 

0.067 

-0.065 

32 

2.141 

0.005 

0.000 

0.005 

11 

1.065 

0.008 

0.147 

-0.138 

28 

1.936 

0.009 

0.005 

0.004 

12 

1.116 

0.060 

0.264 

-0.204 

33 

2.192 

0.004 

0.000 

0.004 

13 

1.167 

0.289 

0.408 

-0.119 

34 

2.243 

0.003 

0.000 

0.003 

14 

1.218 

0.599 

0.557 

0.042 

35 

2.295 

0.002 

0.000 

0.002 

15 

1.270 

0.811 

0.686 

0.124 

36 

2.346 

0.002 

0.000 

0.002 

16 

1.321 

0.913 

0.768 

0.145 

37 

2.397 

0.001 

0.000 

0.001 

17 

1.372 

0.949 

0.788 

0.161 

38 

2.446 

0.001 

0.000 

0.001 

18 

1.423 

0.923 

0.741 

0.182 

39 

2.500 

0.001 

0.000 

0.001 

19 

1.475 

0.741 

0.637 

0.104 

40 

2.551 

0.001 

0.000 

0.001 

20 

1.526 

0.470 

0.502 

-0.032 

27 

1.885 

0.010 

0.010 

0 . 000 

21 

1.577 

0.240 

0.363 

-0.123 

1 

0.552 

0.000 

0.000 

0.000 

22 

1.628 

0.117 

0.242 

-0.124 

2 

0.603 

0.000 

0.000 

0.000 

23 

1 .680 

0.052 

0.147 

-0.095 

3 

0.655 

0.000 

0.000 

0.000 

24 

1.731 

0.026 

0.084 

-0.058 

4 

0.706 

0.000 

0.000 

0.000 

25 

1.782 

0.016 

0.044 

-0.028 

5 

0.757 

0.000 

0.000 

0.000 

26 

1.833 

0.012 

0.022 

-0.010 

6 

0.808 

0.000 

0.000 

0.000 

27 

1 . 885 

0.010 

0.010 

0.000 

7 

0.860 

0.000 

0.001 

-0.001 

28 

1.936 

0.009 

0.005 

0.004 

8 

0.911 

0.000 

0.007 

-0.007 

29 

1.987 

0.008 

0.002 

0.006 

26 

1.833 

0.012 

0.022 

-0.010 

30 

2 . 038 

0.007 

0.001 

0.006 

9 

0.962 

0.001 

0.025 

-0.024 

31 

2.090 

0.006 

0.000 

0.006 

25 

1.782 

0.016 

0.044 

-0.028 
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32  2.141  0.005  0.000  0.005 

33  2.192  0.004  0.000  0.004 

34  2.243  0.003  0.000  0.003 

35  2.295-  0.002  0.000  0.002 

36  2.346  0.002  0.000  0.002 

37  2.397  0.001  0.000  0.001 

38  2.448  0.001  0.000  0.001 

39  2.500  0.001  0.000  0.001 

40  2.551  0.001  0.000  0.001 


END  OF  PROBLEM 


20  1.526  0.470  0.502  -0.032 

24  1.731  0.026  0.084  -0.058 

10  1.013  0.002  0.067  -0.065 

23  1.680  0.052  0.147  -0.095 

13  1.167  0.289  0.408  -0.119 

21  1.577  0.240  0.363  -0.123 

22  1.628  0.117  0.242  -0.124 

11  1.065  0.008  0.147  -0.138 

12  1.116  0.060  0.264  -0.204 
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NON-LINEAR  LEAST  SQUARES  ANALYSIS 

EQUILIBRIUM  TRANSPORT  (MODEL  A) 
THIRD-TYPE  BOUNDARY  CONDITION 
PCE  Run:  1%  Methanol/99%  Water 


INITIAL  VALUES  OF  COEFFICIENTS 


name  initial  value 

PECLET  153.000 

RF  1.000 

PUL  0.337 


OBSERVED  DATA 


OBS.  NO. 
1 
2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 


PORE  VOLUME 
0.5040 
0.5550 
0.6070 
0.6580 
0.7090 
0.7600 
0.8120 
0.8630 
0.9140 
0.9650 
1.0170 
1.0680 
1.1190 
1.1700 
1.2220 
1.2730 
1.3240 
1.3750 
1.4270 
1.4780 
1.5290 
1,5800 
1.6320 
1.6830 
1.7340 
1.7850 
1.8370 
1.8880 
1.9390 
1.9900 
2.0420 
2.0930 
2.1440 
2.1950 
2,2470 
2.2980 
2.3490 
2.4000 


CONCENTRATION 

0.0000 

0.0001 

0.0005 

0.0013 

0.0028 

0.0057 

0.0099 

0.0165 

0.0247 

0.0359 

0.0530 

0.0820 

0.1413 

0.2641 

0.4469 

0.6510 

0.8571 

0.8900 

0.8853 

0.7760 

0.6292 

0.4164 

0.2006 

0.0725 

0.0227 

0.0105 

0.0073 

0.0059 

0.0052 

0.0050 

0.0048 

0.0047 

0.0048 

0.0050 

0.0049 

0.0048 

0.0042 

0.0038 
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39 

40 


2.4520  0.0034 

2.5030  0.0031 


1 


ITERATION 

0 

1 

2 

3 

4 

5 

CORRELATION 


SSQ 

3.5295768 

0.4832896 

0.1273274 

0.1200646 

0.1184520 

0.1184520 

MATRIX 


RF 

1.00000 

1.15290 

1.23162 

1.21627 

1.22041 

1.22041 


1 

1  1.0000 


NON-LINEAR  LEAST  SQUARES  ANALYSIS,  FINAL  RESULTS 


VARIABLE  NAME 
1  RF 


VALUE  S.E.COEFF. 

1.22041  0.0071 


95V  CONFIDENCE  LIMITS 
T-VALUE  LOWER  UPPER 

171.67  1.2060  1.2348 


ORDERED  BY  COMPUTER  INPUT 


PORE  CONCENTRATION  RESI- 

NO  VOLUME  OBS .  FITTED 

DUAL 

PORE 

NO 

1 

0.504 

0.000 

0.000 

0,000 

19 

2 

0.555 

0.000 

0.000 

0.000 

"  17 

3 

0.607 

0.001 

0.000 

0.000 

18 

4 

0.658 

0.001 

0.000 

0.001 

20 

5 

0.709 

0.003 

0.000 

0.003 

21 

6 

0.760 

0.006 

0.000 

0.006 

9 

7 

0.812 

0.010 

0.000 

0.010 

16 

8 

0.863 

0.017 

0.001 

0.015 

10 

9 

0.914 

0.025 

0.005 

0.019 

8 

10 

0.965 

0.036 

0.019 

0.016 

7 

11 

1.017 

0.053 

0.055 

-0.002 

6 

12 

1.068 

0.082 

0.121 

-0.039 

34 

13 

1.119 

0.141 

0.223 

-0.082 

35 

14 

1.170 

0.264 

0.355 

-0.091 

36 

15 

1.222 

0.447 

0.502 

-0.055 

33 

16 

1.273 

0.651 

0.634 

0.017 

37 

17 

1.324 

0.857 

0,732 

0.125 

32 

18 

1.375 

0.890 

0.775 

0.115 

38 

19 

1.427 

0.885 

0.755 

0.131 

39 

20 

1.478 

0.776 

0.676 

0.100 

31 

21 

1.529 

0.629 

0.558 

0.071 

40 

22 

1.580 

0.416 

0.425 

-0.006 

5 

23 

1.632 

0.201 

0.296 

-0.096 

4 

24 

1.683 

0.072 

0.193 

-0.120 

30 

25 

1.734 

0.023 

0.117 

-0.094 

3 

26 

1.785 

0.011 

0.066 

-0.056 

2 

27 

1.837 

0.007 

0.035 

-0.027 

1 

28 

1.888 

0,006 

0.017 

-0.012 

11 

29 

1.939 

0.005 

0.008 

-0.003 

29 

30 

1.990 

0.005 

0.004 

0.001 

22 

31 

2.042 

0.005 

0.002 

0.003 

28 

- ORDERED  BY  RESIDUALS- - 

CONCENTRATION  RES I - 
VOLUME  OBS .  FITTED  DUAL 


1.427 

0.885 

1.324 

0.857 

1.375 

0.890 

1.478 

0.776 

1.529 

0.629 

0.914 

0.025 

1.273 

0.651 

0.965 

0.036 

0.863 

0.017 

0.812 

0.010 

0.760 

0.006 

2.195 

0.005 

2.247 

0.005 

2.298 

0.005 

2.144 

0.005 

2.349 

0.004 

2.093 

0.005 

2.400 

0.004 

2.452 

0.003 

2.042 

0.005 

2.503 

0.003 

0.709 

0.003 

0.656 

0.001 

1.990 

0.005 

0.607 

0.001 

0.555 

0.000 

0.504 

0.000 

1.017 

0.053 

1.939 

0.005 

1.580 

0.416 

1.888 

0.006 

0.755 

0.131 

0.732 

0.125 

0.775 

0.115 

0.676 

0.100 

0.558 

0.071 

0.005 

0.019 

0.634 

0.017 

0.019 

0.016 

0.001 

0.015 

0.000 

0.010 

0.000 

0.006 

0.000 

0.005 

0.000 

0.005 

0.000 

0.005 

0.000 

0.005 

0.000 

0.004 

0.001 

0.004 

0.000 

0.004 

0.000 

0.003 

0.002 

0.003 

0.000 

0.003 

0.000 

0.003 

0.000 

0,001 

0.004 

0.001 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.055 

-0.002 

0.008 

-0.003 

0.425 

-0.008 

0.017 

-0.012 
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32  2.093  0.005  0.001  0.004 

33  2.144  0.005  0.000  0.005 

34  2.195  0.005  0.000  0.005 

35  2.247'  0.005  0.000  0.005 

36  2.298  0.005  0.000  0.005 

37  2.349  0.004  0.000  0.004 

38  2.400  0.004  0.000  0.004 

39  2.452  0.003  0.000  0.003 

40  2.503  0.003  0.000  0.003 


END  OF  PROBLEM 


27 

1.837 

0.007 

0.035 

-0.027 

12 

1.068 

0.082 

0.121 

-0.039 

15 

1.222 

0.447 

0.502 

-0.055 

26 

1.785 

0.011 

0.066 

-0.056 

13 

1.119 

0.141 

0.223 

-0.082 

14 

1.170 

0.264 

0.355 

-0.091 

25 

1.734 

0.023 

0.117 

-0.094 

23 

1.632 

0.201 

0.296 

-0.096 

24 

1.683 

0.072 

0.193 

-0.120 

100 


NON-LINEAR  LEAST  SQUARES  ANALYSIS 

EQUILIBRIim  TRANSPORT  (MODEL  A) 
THIRD-TYPE  BOUNDARY  CONDITION 
PCE  Run:  5%  Methanol/9S%  Water 


INITIAL 

VALUES  OF 

COEFFICIENTS 

NO 

NAME 

INITIAL  VALUE 

1 

PECLET 

305.000 

2 

RF 

1.000 

3 

PUL 

0.309 

OBSERVED 

DATA 

OBS.  NO. 

PORE  VOLUME 

CONCENTRATION 

1 

0.6060 

0.0003 

2 

0.6580 

0.0007 

3 

0.7090 

0.0019 

4 

0.7600 

0.0041 

5 

0.8110 

0.0077 

6 

0.8630 

0.0133 

7 

0.9140 

0.0218 

8 

0.9650 

0.0338 

9 

1.0160 

0.0499 

10 

1.0680 

0.0785 

11 

1.1190 

0.1387 

12 

1.1700 

0.2247 

13 

1.2210 

0.3449 

14 

1.2730 

0.5036 

15 

1,3240 

0.7017 

16 

1.3750 

0.8171 

17 

1.4260 

0.8316 

18 

1.4780 

0.7858 

19 

1.5290 

0.6491 

20 

1.5800 

0.3957 

21 

1.6310 

0.1656 

22 

1.6830 

0.0522 

23 

1.7340 

0.0200 

24 

1,7850 

0.0109 

25 

1.8360 

0.0070 

26 

1.8880 

0.0055 

27 

1.9390 

0.0047 

28 

1.9900 

0.0041 

29 

2.0410 

0.0038 

30 

2.0930 

0.0036 

31 

2.1440 

0.0034 

32 

2.1950 

0.0032 

33 

2.2460 

0.0031 

34 

2,2980 

0.0029 

35  . 

2.3490 

0.0028 

36 

2.4000 

0.0027 

37 

2.4510 

0.0026 

38 

2.5030 

0.0027 
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39 

40 


2.5540  0.0025 

2,6050  0.0023 


1 


ITERATION 

0 

1 

2 

3 

4 

5 

CORRELATION 


SSQ 

4.4736667 

2.1920920 

0.3298729 

0.0400964 

0.0381951 

0.0381951 

MATRIX 


RF 

1.00000 

1.09357 

1.19663 

1.24727 

1.25170 

1.25171 


1 

1  1.0000 


NON-LINEAR  LEAST  SQUARES  ANALYSIS,  FINAL  RESULTS 


VARIABLE  NAME 
1  RF 


VALUE  S . E . COEFF . 

1.25171  0.0032 


95%  CONFIDENCE  LIMITS 
T- VALUE  LOWER  UPPER 

394.36  1.2453  1.2581 


- ORDERED  BY  COMPUTER  INPUT . 

PORE  CONCENTRATION  RES I-  PORE 


NO 

VOLUME 

OBS. 

1 

0.606 

0.000 

2 

0.658 

0.001 

3 

0.709 

0.002 

4 

0.760 

0.004 

5 

0.811 

0.008 

6 

0.863 

0.013 

7 

0.914 

0.022 

8 

0.965 

0.034 

9 

1.016 

0.050 

10 

1.068 

0.079 

11 

1.119 

0.139 

12 

1.170 

0.225 

13 

1.221 

0.345 

14 

1.273 

0.504 

15 

1.324 

0.702 

16 

1.375 

0.817 

17 

1.426 

0.832 

18 

1.478 

0.786 

19 

1.529 

0.649 

20 

1.580 

0.396 

21 

1.631 

0.166 

22 

1.683 

0.052 

23 

1.734 

0.020 

24 

1.785 

0.011 

25 

1.836 

0.007 

26 

1.888 

0.005 

27 

1.939 

0.005 

28 

1.990 

0.004 

29 

2.041 

0.004 

30 

2.093 

0.004 

31 

2.144 

0,003 

FITTED 

DUAL 

NO 

0.000 

0.000 

11 

0.000 

0.001 

10 

0.000 

0.002 

9 

0.000 

0.004 

8 

0.000 

0.008 

19 

0.000 

0.013 

12 

0.000 

0.022 

7 

0.001 

0.033 

6 

0.005 

0.045 

5 

0.025 

0.054 

18 

0.083 

0.056 

27 

0.202 

0.023 

4 

0.379 

-0.034 

28 

0.582 

-0.078 

29 

0.752 

-0.050 

30 

0.854 

-0.037 

26 

0.867 

-0.036 

31 

0.781 

0.004 

32 

0.618 

0.031 

33 

0.423 

-0.027 

34 

0.249 

-0.084 

35 

0.124 

-0.072 

38 

0.054 

-0.034 

36 

0.021 

-0.010 

37 

0.007 

0.000 

39 

0.002 

0.004 

40 

0.001 

0.004 

3 

0.000 

0.004 

2 

0.000 

0.004 

1 

0.000 

0.004 

25 

0.000 

0.003 

24 

. ORDERED  BY  RESIDUALS- - 

CONCENTRATION  RES I - 
VOLUME  OBS .  FITTED  DUAL 


1.119 

0.139 

1.068 

0.079 

1.016 

0.050 

0.965 

0.034 

1.529 

0.649 

1.170 

0.225 

0.914 

0.022 

0.863 

0.013 

0.811 

0.008 

1.478 

0.786 

1.939 

0.005 

0.760 

0.004 

1.990 

0.004 

2.041 

0.004 

2.093 

0.004 

1.888 

0.005 

2.144 

0.003 

2.195 

0.003 

2.246 

0.003 

2.298 

0.003 

2.349 

0.003 

2.503 

0.003 

2.400 

0.003 

2.451 

0.003 

2.554 

0.003 

2.605 

0.002 

0.709 

0.002 

0.658 

0.001 

0.606 

0.000 

1.836 

0.007 

1.785 

0.011 

0.083 

0.056 

0.025 

0.054 

0.005 

0.045 

0.001 

0.033 

0.618 

0.031 

0.202 

0.023 

0.000 

0.022 

0.000 

0.013 

0.000 

0.008 

0.781 

0.004 

0,001 

0.004 

0.000 

0.004 

0.000 

0.004 

0.000 

0.004 

0.000 

0.004 

0.002 

0.004 

0.000 

0.003 

0.000 

0.003 

0.000 

0.003 

0,000 

0.003 

0.000 

0.003 

0.000 

0.003 

0.000 

0.003 

0,000 

0.003 

0.000 

0.002 

0.000 

0.002 

0.000 

0.002 

0.000 

0.001 

0.000 

0.000 

0.007 

0.000 

0,021 

-0.010 
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32  2.195  0.003  0.000  0.003 

33  2.246  0.003  0.000  0.003 

34  2.298  0.003  0.000  0.003 

35  2.349-  0.003  0.000  0.003 

36  2.400  0.003  0.000  0,003 

37  2.451  0.003  0.000  0.003 

38  2.503  0.003  0.000  0.003 

39  2.554  0.003  0.000  0.002 

40  2.605  0.002  0,000  0.002 


END  OF  PROBLEM 


20  1.580  0.396  0.423  -0.027 

23  1.734  0.020  0.054  -0.034 

13  1.221  0.345  0.379  -0.034 

17  1.426  0.832  0.867  -0.036 

16  1.375  0.817  0,854  -0.037 

15  1.324  0.702  0.752  -0,050 

22  1.683  0.052  0.124  -0.072 

14  1.273  0.504  0.582  -0.078 

21  1.631  0.166  0.249  -0.084 
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NON-LINEAR  LEAST  SQUARES  ANALYSIS 


EQUILIBRIUM  TRANSPORT  (MODEL  A) 
THIRD-TYPE  BOUNDARY  CONDITION 
PCE  Run:  lOV  Methanol/90%  Water 


INITIAL  VALUES  OF  COEFFICIENTS 


NO 

NAME 

INITIAL  VALUE 

1 

PECLET 

276.000 

2 

RF 

1.000 

3 

PUL 

0.320 

OBSERVED 

DATA 

DBS .  ■  NO . 

PORE  VOLUME 

CONCENTRATION 

1 

0.4500 

0.0002 

2 

0.5010 

0.0007 

3 

0.5530 

0.0018 

4 

0.6040 

0.0044 

5 

0.6550 

0.0077 

6 

0.7060 

0.0122 

7 

0.7580 

0.0181 

6 

0.8090 

0.0242 

9 

0.8600 

0.0328 

10 

0.9110 

0.0406 

11 

0.9630 

0.0507 

12 

1.0140 

0.0659 

13 

1.0650 

0.0934 

14 

1.1160 

0.1523 

15 

1.1680 

0.2284 

16 

1.2190 

0.3264 

17 

1.2700 

0.5001 

18 

1.3210 

0.6810 

19 

1.3730 

0.7971 

20 

1.4240 

0.8124 

21 

1.4750 

0.7504 

22 

1.5260 

0.6556 

23 

1.5780 

0.4617 

24 

1.6290 

0.2220 

25 

1.6800 

0.0719 

26 

1.7310 

0.0266 

27 

1.7830 

0.0139 

28 

1.8340 

0.0093 

29 

1.8850 

0.0068 

30 

1.9360 

0.0054 

31 

1.9880 

0.0045 

32 

2.0390 

0.0038 

33 

2.0900 

0.0033 

34 

2.1410 

0.0029 

35 

2.1930 

0.0026 

36 

2.2440 

0.0024 

37 

2.2950 

0.0022 

38 

2.3460 

0.0021 
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39 

40 


2,3980  0.0020 

2.4490  0.0019 


1 


ITERATION 

0 

1 

2 

3 

4 

5 

CORRELATION 


SSQ 

4.1831929 

1.9238347 

0.2849047 

0.0550446 

0.0528172 

0.0528172 

MATRIX 


RF 

1.00000 

1.09672 

1.19765 

1.24384 

1.24887 

1.24887 


1 

1  1.0000 


NON-LINEAR  LEAST  SQUARES  ANALYSIS,  FINAL  RESULTS 


VARIABLE  NAME 
1  RF 


VALUE  S . E . COEFF . 

1.24887  0.0038 


95%  CONFIDENCE  LIMITS 
T- VALUE  LOWER  UPPER 

325.63  1.2411  1.2566 


. ORDERED  BY  COMPUTER 

PORE  CONCENTRATION  RES I- 

INPUT-- 

PORE 

NO 

VOLUME 

OBS. 

FITTED 

DUAL 

NO 

1 

0.450 

0.000 

0.000 

0.000 

13 

2 

0.501 

0.001 

0.000 

0.001 

14 

3 

0.553 

0.002 

0.000 

0.002 

12 

4 

0.604 

0.004 

0.000 

0.004 

11 

5 

0.655 

0.008 

0.000 

0.008 

10 

6 

0.706 

0.012 

0.000 

0.012 

9 

7 

0.758 

0.018 

0.000 

0.018 

8 

8 

0.809 

0.024 

0.000 

0.024 

7 

9 

0.860 

0.033 

0.000 

0.033 

15 

10 

0.911 

0.041 

0.000 

0.041 

6 

11 

0.963 

0.051 

0.001 

0.050 

5 

12 

1.014 

0.066 

0.007 

0.059 

22 

13 

1.065 

0.093 

0.030 

0.063 

4 

14 

1.116 

0.152 

0.093 

0.060 

30 

15 

1.168 

0.228 

0.215 

0.013 

31 

16 

1.219 

0.326 

0.388 

-0.061 

32 

17 

1.270 

0.500 

0.578 

-0.077 

33 

18 

1.321 

0.681 

0.741 

-0.060 

29 

19 

1.373 

0.797 

0.845 

-0.048 

34 

20 

1.424 

0.812 

0.866 

-0.053 

35 

21 

1.475 

0.750 

0.796 

-0.046 

36 

22 

1.526 

0.656 

0.650 

0.005 

37 

23 

1.578 

0.462 

0.463 

-0.001 

38 

24 

1.629 

0.222 

0.289 

-0.067 

39 

25 

1.680 

0.072 

0.158 

-0.086 

40 

26 

1.731 

0.027 

0.075 

-0.049 

3 

27 

1.783 

0.014 

0.031 

-0.017 

2 

28 

1.834 

0*-009 

0.012 

-0.002 

1 

29 

1.885 

0.007 

0.004 

0.003 

23 

30 

1.936 

0.005 

0.001 

0.004 

28 

31 

1.988 

0.004 

0.000 

0.004 

27 

- ORDERED  BY  RESIDUALS -- 

CONCENTRATION  RES I - 

VOLUME  OBS.  FITTED  DUAL 


1.065 

0.093 

1.116 

0.152 

1.014 

0.066 

0.963 

0.051 

0.911 

0.041 

0.860 

0.033 

0.809 

0.024 

0.758 

0.018 

1.168 

0.228 

0.706 

0.012 

0.655 

0.008 

1.526 

0.656 

0.604 

0.004 

1.936 

0.005 

1.988 

0.004 

2.039 

0.004 

2.090 

0.003 

1.885 

0.007 

2.141 

0.003 

2.193 

0.003 

2.244 

0.002 

2.295 

0.002 

2.346 

0.002 

2.398 

0.002 

2.449 

0.002 

0.553 

0.002 

0.501 

0.001 

0.450 

0.000 

1.578 

0.462 

1.834 

0.009 

1.783 

0.014 

0.030 

0.063 

0.093 

0.060 

0.007 

0.059 

0.001 

0.050 

0.000 

0.041 

0.000 

0.033 

0.000 

0.024 

0.000 

0.018 

0.215 

0.013 

0.000 

0.012 

0.000 

0.008 

0.650 

0.005 

0.000 

0,004 

0.001 

0.004 

0.000 

0.004 

0.000 

0.004 

0.000 

0.003 

0.004 

0.003 

0.000 

0.003 

0.000 

0.003 

0.000 

0.002 

0.000 

0.002 

0.000 

0.002 

0.000 

0.002 

0.000 

0.002 

0.000 

0.002 

0.000 

0.001 

0.000 

0.000 

0.463 

-0.001 

0.012 

-0.002 

0.031 

-0.017 
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32 

2,039 

0.004 

33 

2.090 

0.003 

34 

2.141 

0.003 

35 

2.193  ■ 

0.003 

36 

2.244 

0.002 

37 

2.295 

0.002 

30 

2.346 

0.002 

39 

2.398 

0.002 

40 

2.449 

0.002 

0.000 

0.004 

21 

0.000 

0.003 

19 

0.000 

0.003 

26 

0.000 

0.003 

20 

0.000 

0.002 

18 

0.000 

0.002 

16 

0.000 

0.002 

24 

0.000 

0.002 

17 

0.000 

0.002 

25 

1.475 

0.750 

0.796 

1.373 

0.797 

0.845 

1.731 

0.027 

0.075 

1.424 

0.812 

0.866 

1.321 

0.681 

0.741 

1.219 

0.326 

0.386 

1.629 

0.222 

0.289 

1.270 

0.500 

0.578 

1.680 

0.072 

0.158 

END  OF  PROBLEM 


0,046 

0.048 

0.049 

0.053 

0.060 

0,061 

0.067 

0.077 

0.086 
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1 


NON-LINEAR  LEAST  SQUARES  ANALYSIS 


EQUILIBRIUM  TRANSPORT  (MODEL  A) 
THIRD-TYPE  BOUNDARY  CONDITION 
PCE  Run:  20%  Methanol/80%  Water 


INITIAL  VALUES  OF  COEFFICIENTS 


NO  NAME  INITIAL  VALUE 

1  PECLET  328.000 

2  RF  1.000 

3  PUL  0.313 


OBSERVED  DATA 


OBS.  NO. 
1 
2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 
35. 

36 

37 

38 


PORE  VOLUME 
0.4830 
0.5350 
0.5860 
0.6370 
0.6880 
0.7400 
0.7910 
0.8420 
0.8930 
0.9450 
0.9960 
1.0470 
1.0980 
1.1500 
1.2010 
1.2520 
1.3030 
1.3550 
1.4060 
1.4570 
1.5080 
1.5600 
1.6110 
1.6620 
1.7130 
1.7650 
1.8160 
1.8670 
1.9180 
1.9700 
2.0210 
2.0720 
2.1230 
2.1750 
2.2260 
2.2770 
2.3280 
2.3800 


CONCENTRATION 

0.0125 

0.0172 

0.0221 

0.0276 

0.0330 

0.0374 

0.0444 

0.0516 

0.0607 

0.0735 

0.0930 

0.1491 

0.2285 

0.3114 

0.3575 

0.4748 

0.6591 

0.8305 

0.8706 

0.7595 

0.6463 

0.4454 

0.2223 

0.0708 

0.0333 

0.0120 

0.0140 

0.0111 

0.0092 

0.0075 

0.0072 

0.0056 

0.0052 

0.0047 

0.0046 

0.0044 

0.0043 

0.0044 
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39 

2.4310 

0.0043 

40 

2.4820 

0.0040 

ITERATION 

SSQ 

RF 

0 

3.7921204 

1.00000 

1 

1.8465650 

1.08619 

2 

0.4169640 

1.17352 

3 

0.1260098 

1.22018 

4 

0.1150433 

1.23034 

5 

0.1150326 

1.23175 

6 

0.1150326 

1.23175 

CORRELATION 

MATRIX 

1 

1  1.0000 

NON-LINEAR  LEAST  SQUARES  ANALYSIS,  FINAL  RESULTS 


VARIABLE  NAME 
1  RF 


95%  CONFIDENCE  LIMITS 

VALUE  S.E.COEFF.  T- VALUE  LOWER  UPPER 

1.23175  0.0053  233.02  1.2211  1.2424 


. ORDERED  BY  COMPUTER  INPUT 

PORE  CONCENTRATION  RESI- 


NO 

VOLUME 

OBS. 

FITTED 

DUAL 

1 

0.483 

0.013 

0.000 

0.013 

2 

0.535 

0.017 

0.000 

0.017 

3 

0.586 

0.022 

0.000 

0.022 

4 

0.637 

0.028 

0.000 

0.028 

5 

0.688 

0.033 

0.000 

0.033 

6 

0.740 

0.037 

0.000 

0.037 

7 

0.791 

0.044 

0.000 

0.044 

8 

0.842 

0.052 

0.000 

0.052 

9 

0.893 

0.061 

0.000 

0.061 

10 

0.945 

0.073 

0.000 

0,073 

11 

0.996 

0.093 

0.003 

0.090 

12 

1.047 

0.149 

0.018 

0.131 

13 

1.098 

0.229 

0.070 

0.158 

14 

1.150 

0.311 

0.189 

0.122 

15 

1.201 

0.357 

0.373 

-0.015 

16 

1.252 

0.475 

0.583 

-0.108 

17 

1.303 

0.659 

0.762 

-0.103 

18 

1.355 

0.831 

0.874 

-0.043 

19 

1.406 

0.871 

0.893 

-0.022 

20 

1.457 

0.759 

0.813 

-0.053 

21 

1.508 

0.646 

0.647 

0.000 

22 

1.560 

0.445 

0.436 

0.009 

23 

1.611 

0.222 

0.251 

-0.028 

24 

1.662 

0.071 

0.122 

-0.051 

25 

1.713 

0.033 

0.050 

-0.017 

26 

1.765 

0.012 

0.017 

-0,005 

27 

1.816 

0.014 

0.005 

0.009 

28 

1.867 

0,011 

0.001 

0.010 

29 

1.918 

0.009 

0.000 

0.009 

30 

1.970 

0.007 

0.000 

0.007 

ORDERED  BY  RESIDUALS 


PORE 

CONCENTRATION 

RES  I 

- 

NO 

VOLUME 

OBS. 

FITTED 

DUAL 

13 

1.098 

0.229 

0.070 

0.158 

12 

1.047 

0.149 

0.018 

0.131 

14 

1.150 

0.311 

0.189 

0.122 

11 

0.996 

0.093 

0.003 

0.090 

10 

0.945 

0.073 

0.000 

0.073 

9 

0.893 

0.061 

0.000 

0.061 

8 

0.842 

0.052 

0.000 

0.052 

7 

0.791 

0.044 

0.000 

0.044 

6 

0.740 

0.037 

0.000 

0.037 

5 

0.688 

0.033 

0.000 

0.033 

4 

0.637 

0.028 

0.000 

0.028 

3 

0.586 

0.022 

0.000 

0.022 

2 

0.535 

0.017 

0.000 

0,017 

1 

0.483 

0,013 

0.000 

0.013 

28 

1.867 

0.011 

0.001 

0.010 

22 

1.560 

0.445 

0.436 

0.009 

29 

1.918 

0.009 

0.000 

0.009 

27 

1.816 

0.014 

0.005 

0.009 

30 

1.970 

0.007 

0.000 

0.007 

31 

2.021 

0.007 

0.000 

0.007 

32 

2.072 

0.006 

0.000 

0.006 

33 

2.123 

0.005 

0.000 

0.005 

34 

2.175 

0.005 

0.000 

0.005 

35 

2.226 

0.005 

0.000 

0.005 

36 

2.380 

0.004 

0.000 

0.004 

36 

2.277 

0.004 

0.000 

0.004 

39 

2.431 

0.004 

0.000 

0.004 

37 

2.328 

0.004 

0.000 

0.004 

40 

2.482 

0.004 

0.000 

0.004 

21 

1.508 

0.646 

0.647 

0.000 
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31 

2.021 

0.007 

0.000 

0.007 

32 

2.072 

0.006 

0.000 

0.006 

33 

2.123 

0.005 

0.000 

0.005 

34 

2.175- 

0.005 

0.000 

0.005 

35 

2.226 

0.005 

0.000 

0.005 

36 

2.277 

0.004 

0.000 

0.004 

37 

2.328 

0.004 

0.000 

0.004 

38 

2.380 

0.004 

0.000 

0.004 

39 

2.431 

0.004 

0.000 

0.004 

40 

2.482 

0.004 

0.000 

0.004 

26 

1.765 

0.012 

0.017 

-0.005 

15 

1.201 

0.357 

0.373 

-0.015 

25 

1.713 

0.033 

0.050 

-0.017 

19 

1.406 

0.871 

0.893 

-0.022 

23 

1.611 

0.222 

0.251 

-0.028 

18 

1.355 

0.831 

0.874 

-0.043 

24 

1.662 

0.071 

0.122 

-0.051 

20 

1.457 

0.759 

0.813 

-0.053 

17 

1.303 

0.659 

0.762 

-0.103 

16 

1.252 

0.475 

0.563 

-0.108 

END  OF  PROBLEM 


1 


NON-LINEAR  LEAST  SQUARES  ANALYSIS 

EQUILIBRIUM  TRANSPORT  (MODEL  A) 
THIRD-TYPE  BOUNDARY  CONDITION 
PCE  Run:  50V  Methanol/50%  Water 


INITIAL  VALUES  OF  COEFFICIENTS 


10  NAME  INITIAL  VALUE 

1  PECLET  221.000 

2  RF  1.000 

3  PUL  0.326 


OBSERVED  DATA 


OBS.  NO. 
1 
2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 
35- 

36 

37 

38 


PORE  VOLUME 
0.4330 
0.4850 
0.5360 
0.5870 
0.6380 
0.6900 
0.7410 
0.7920 
0.8430 
0.8950 
0.9460 
0.9970 
1.0480 
1.1000 
1.1510 
1.2020 
1.2530 
1.3050 
1.3560 
1.4070 
1.4580 
1.5100 
1.5610 
1.6120 
1.6630 
1.7150 
1.7660 
1.8170 
1.8680 
1.9200 
1.9710 
2.0220 
2.0730 
2.1250 
2.1760 
2.2270 
2.2780 
2.3300 


CONCENTRATION 

0.0084 

0.0136 

0.0191 

0.0237 

0.0296 

0.0352 

0.0384 

0.0453 

0.0511 

0.0757 

0.0688 

0.0902 

0.1298 

0.2016 

0.2740 

0.3354 

0.4179 

0.5865 

0.8007 

0.9018 

0.8672 

0.7320 

0.4621 

0.2397 

0.1438 

0.0839 

0.0619 

0.0458 

0.0300 

0.0191 

0.0122 

0.0083 

0.0057 

0.0048 

0.0039 

0.0034 

0.0030 

0.0026 


no 


3  9 

2.3810 

0.0024 

40 

2.4320 

0.0022 

ITERATION 

SSQ 

RF 

0 

3.8082657 

1.00000 

1 

1.5077115 

1.10596 

2 

0.1944039 

1.20882 

3 

0.1126366 

1.23999 

4 

0.1126193 

1.24035 

CORRELATION 

MATRIX 

1 

1  1.0000 


non-linear  least  squares  analysis,  final  resets 


VARIABLE  NAME 
1  RF 


VALUE  S . E . COEFF . 

1.24035  0.0061 


95%  CONFIDENCE  LIMITS 

t-value  lower  upper 

204.98  1.2281  1.2526 


. ordered  by  computer 

PORE  CONCENTRATION  RESI- 
NO  VOLUME  OBS .  FITTED 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 


433 

485 

536 

587 

638 

690 

741 

792 

843 

895 

.946 

,997 

,048 

.100 

.151 

.202 

.253 

.305 

.356 

.407 

.458 

.510 

.561 

.612 

.663 

.715 

.766 

.817 

.868 

.920 

.971 

.022 


008 

014 

019 

024 

030 

035 

038 

045 

051 

.076 

,069 

,090 

,130 

.202 

.274 

.335 

.418 

.587 

.801 

.902 

.867 

.732 

.462 

.240 

.144 

.084 

.062 

046 


0, 

0, 

0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 

0.030 

0.019 


.012 

.008 


000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

002 

oil 

038 

,103 

,215 

,370 

.541 

.698 

.801 

.835 

.788 

.669 

.511 

.349 

.214 

.116 

.058 

.026 

.011 

.004 

.001 

.000 


INPUT- 

DUAL 
0.008 
0.014 
0.019 
0.024 
0.030 
0.035 
0.038 
0,045 
0.051 
0.075 
0.067 
0.080 
092 
099 
059 
035 
124 
111 
0.001 
0.067 
0.079 
0.063 
>0.049 
-0.109 
070 
032 
004 
020 
019 
015 
Oil 
008 


-0. 

-0. 

0. 

0. 

0. 

0. 

0. 

0. 


PORE 

NO 

14 
13 
12 
21 
10 
20 
11 
22 

15 
9 
8 
7 
6 
5 
4 

28 

29 
3 

30 
2 

31 
1 

32 

33 

34 
27 

35 

36 

37 

38 

39 

40 


■ORDERED  BY  RESIDUALS 


CONCENTRATION 


VOLUME 

OBS. 

1.100 

0.202 

1.048 

0.130 

0.997 

0.090 

1.458 

0.867 

0.895 

0.076 

1.407 

0.902 

0.946 

0.069 

1.510 

0.732 

1.151 

0.274 

0.843 

0.051 

0.792 

0.045 

0.741 

0.038 

0.690 

0.035 

0.638 

0.030 

0.587 

0.024 

1.817 

0.046 

1.868 

0.030 

0.536 

0.019 

1.920 

0.019 

0.485 

0.014 

1.971 

0.012 

0.433 

0.008 

2.022 

0.008 

2.073 

0.006 

2,125 

0.005 

1.766 

0.062 

2.176 

0.004 

2.227 

0.003 

2.278 

0.003 

2.330 

0.003 

2.381 

0.002 

2.432 

0.002 

RESI- 

*ed  dual 


0.103 

0.099 

0.038 

0.092 

0.011 

0.080 

0.788 

0.079 

0.000 

0.075 

0.835 

0.067 

0.002 

0.067 

0.669 

0.063 

0.215 

0,059 

0.000 

0.051 

0.000 

0.045 

0.000 

0.038 

0.000 

0.035 

0.000 

0.030 

0.000 

0.024 

0.026 

0.020 

0.011 

0.019 

0.000 

0.019 

0.004 

0.015 

0.000 

0.014 

0.001 

0.011 

0.000 

0.008 

0.000 

0.008 

0.000 

0.006 

0.000 

0.005 

0.058 

0.004 

0.000 

0.004 

0.000 

0.003 

0.000 

0,003 

0.000 

0.003 

0.000 

0.002 

0.000 

0.002 

111 


33 

2.073 

0.006 

34 

2.125 

0.005 

35 

2.176 

0.004 

36 

2.221' 

0.003 

37 

2.278 

0.003 

38 

2.330 

0.003 

39 

2.381 

0.002 

40 

2.432 

0.002 

0.000 

0.006 

19 

0.000 

0.005 

26 

0.000 

0.004 

16 

0.000 

0.003 

23 

0.000 

0.003 

25 

0.000 

0.003 

24 

0.000 

0.002 

18 

0.000 

0.002 

17 

1.356 

0.801 

0.801 

1.715 

0.084 

0.116 

1.202 

0.335 

0.370 

1.561 

0.462 

0.511 

1.663 

0.144 

0.214 

1.612 

0.240 

0.349 

1.305 

0.587 

0.698 

1.253 

0.418 

0.541 

END  OF  PROBLEM 


0.001 

0.032 

0.035 

0.049 

0.070 

0.109 

0.111 

0.124 
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